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INTRODUCTION 



The corner-stone of modern cosmology is that, at least on large scales, the visible universe 
seems to be the same in all directions around us and around all points, i.e. the Universe is 
almost homogeneous and isotropic. This is borne out by a variety of observations, particu- 
lary observations of cosmic microwave background (CMB); this radiation has been traveling 
to us for about 14000 million years (see Fig. II. ip . supporting the conclusion that the Uni- 
verse at sufficiently large distances is nearly the same. On the other hand, it is apparent 
that nearby regions of the observable Universe are at present highly inhomogeneous, with 
material clumped into stars, galaxies and galaxy clusters. It is believed that these struc- 
tures have formed over the time via gravitational attraction, from a distribution that was 
more homogeneous in the past. 

The large-scale behavior of the Universe can be described by assuming a homogeneous 
background. On this background, we can superimpose the short scale irregularities. For 
much of the evolution of the observable Universe, these irregularities can be considered to 
be small perturbations on the evolution of the background (unperturbed) Universe. The 
metric of unperturbed Universe is called the Friedman-Leimatre-Roberson- Walker metric, 
and its line element can be to written as: 

ds"^ = -df + a\t) [dr"^ + r'^{de'^ + sin (f^dct?) , (1.1) 

where a{t) is the scale factor and r, 9, are the spherical comoving coordinate^] 

The model described by the above metric is known as the standard cosmological model 
(known also as Big-Bang cosmological model) [661 11711 11721 11731 1216] and is the successful 



particle in this metric have fixed-coordiantes. 



2 



framework that describes the observed properties of the Universe: homogeneity and isotropy 
at large scales, Hubble expansion, almost 14 billion years of evolution in agreement with 
globular clusters and radioactive isotopes dating, cosmic microwave background radiation 
(CMB) confirmed by Penzias and Wilson's discovery in 1965 [HI 1164] . and the relative 
abundances of light elements [51 [TU |[^ 1^116111215^1217] in full agreement with observation. 




Figure 1.1: A representation of the evolution of the universe over 13.7 billion years. The far 
left depicts the earliest moment we can now probe, when a period of " inflation " produced 
a burst of exponential growth in the universe. (Size is depicted by the vertical extent of 
the grid in this graphic.) For the next several billion years, the expansion of the universe 
gradually slowed down as the matter in the universe pulled on itself via gravity. More 
recently, the expansion has begun to speed up again as the repulsive effects of dark energy 
have come to dominate the expansion of the universe. The afterglow light seen by WMAP 
was emitted about 380,000 years after inflation and has traversed the universe largely 
unimpeded since then. The conditions of earlier times are imprinted on this light; it also 
forms a backlight for later developments of the universe (Courtesy of the NASA/WMAP 
Science Team |159] ). 

The introduction of a period of exponential expansion (called inflationary) [H [821 1131] , prior 
to the Big-Bang, brought an elegant solution to the horizon, flatness, and unwanted relics 
problems that were present in the original standard cosmological model [8ll82 [[TT7lll31[[T69] . 
In spite of its success at solving the above mentioned problems, the inflationary period be- 
came perhaps more important because of its ability to stretch the quantum fluctuations 
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of the fields living in the FRW spacetime [HI ESI EH [ml [TMl [TSEl [1691 [2Ql], making 
them classical [3 [361 [781 [H [1331 EM EM [M [ISDl MQ\ and almost constant soon 
after horizon exit. They correspond to small inhomogeneities in the energy density and are 
responsible, via gravitational attraction, of the large-scale structure seen today in the Uni- 
verse. If this scenario turned to be correct, the energy density inhomogeneities should have 
left their trace in the CMB released at the time of recombination. Indeed, the Cosmic Back- 
ground Explorer (COBE) in 1992 [158] found and mapped small anisotropics in the CMB 
temperature of the order of 1 part in 10^ (with average temperature To = 2.725 ± 0.002 K 
[23]), on scales of order thousands of Megaparsecs. With 30 times better angular resolution 
and sensitivity than COBE, the Wilkinson Microwave Anisotropy Probe (WMAP) [159] 
confirmed this picture (see Fig. II. 2p . measuring in turn the cosmological parameters with 
a 1% order precision [119] on scales of order tens of Megaparsecs. The PLANCK satellite 
[591 1206] , launched in may 2009, will be able to refine these observations (see Fig. 11.31 
and 11.41) . With 10 times better angular resolution and sensitivity than WMAP, PLANCK 
promises to determine the temperature anisotropics with a resolution of the order of 1 part 
in 10®, and the cosmological parameters with a 0.1% order precision. 




Figure 1.2: CMB temperature anisotropics as seen by the WMAP satellite (five years 
resuls) [119] . The oval shape is a projection to display the whole sky. The temperature 
anisotropics are found to be of the order of 1 part in 10^. The background temperature 
is To = 2.725 ± 0.002 K; regions at that temperature are in very light blue. The hottest 
regions (in red) correspond to AT ~ 200yuK. The coldest regions (in very dark blue) 
correspond to AT ~ -200/iK (Courtesy of the NASA/WMAP Science Team [T59] ). 

The anisotropics in the CMB temperatur^ ST/Tq are directly related to the perturbation 



^From now on, and unless otherwise stated, the perturbation 5y in any quantity y will be regarded as 
first-order in cosmological perturbation theory. Unperturbed quantities will be denoted by a subscript 



Figure 1.3: Simulation of the CMB temperature anisotropies as seen by the PLANCK 
sateUite. PLANCK will provide a map of the CMB field at all angular resolutions greater 
than 10 arcminutes and with a temperature resolution of the order of 1 part in 10^ (ten 
times better than WMAP) (Courtesy of ESA's PLANCK mission j59]). 

in the spatial curvature C (Sachs- Wolfe effect), whose primarily origin is the stretched 
quantum fluctuations of one or several scalar fields 0j that fill the Universe during inflation 



The quantity ( is related to the perturbation in the intrinsic curvature of space-time slices 
with uniform energy density |149] : 



where ijj is the first order scalar perturbation in the spacial metric. 

Astronomers work with the observable quantity 6T/T and theoretical cosmologists work 
with (. Therefore, we may study the statistical porperties of the observed 6T/T through the 
spectral functions associated with the primordial curvature perturbation whose properties 
are in general model dependent. Knowing the statistical descriptors of ( for some particular 



[TiOllTTGP : 




(1.2) 



(1.3) 



a 



unless otherwise stated. 

^In this and the following expressions the subscripts k stand for the Fourier modes with comoving 
wavenumber k. 
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Figure 1.4: A map of the area of the sky mapped by PLANCK during the first hght survey. 
The colours indicate the magnitude of the deviations of the temperature of the Cosmic 
Microwave Background from its average value (red is hotter and blue is colder). (Courtesy 
of ESA's PLANCK mission 

and well motivated cosmological model proposed for the origin of large scale strucutre, we 
can reject the model or keep it, because some of the statistical descriptors for 6T/T are 
known with good acuracy or at least have an upper bound |119] . 

The statistical properties of the CMB temperature anisotropies can be then described 
in terms of the spectral functions, like the spectrum, bispectrum, trispectrum, etc., of 
the primordial curvature perturbation (. This spectral functions are given in terms of 
other quantities, which have an observational value or an uppper bound. For example, the 
spectrum is parametrized in terms of an amplitude P^^^^, a spectral index and the level 
of statistical anisotropy (7^; the bispectrum 5^ and trispectrum are parametrized in terms 
of products of the spectrum and the quantities /nl, and tnl and (?nl, respectively. As we 
will see in the next chapter, the statistical descriptors /nl, tnl and (^nl are usually called 
levels of non-gaussianity, because non zero values for these quantities imply non-gaussianity 
in the primordial curvature perturbation ( as well in the constrast in the temperature of the 
CMB radiation 6T/T. The non-gaussian characteristics in the CMB are actually present in 
the observation |119] as we will see in more detail in Section 12.41 The status of observation 
can be summarized as foUowfl the spectral amplitude V^;^^^ = (4.957 ± 0.094) X 10 ^ [35] , 
the spectral index = 0.960 ± 0.014 at 2a |119] . the level of non-gaussianity /nl in the 
bispectrum is in the range —9 < f^i < 111 at 2a |119] : and there is no observational 



^We are using values according of the five year of data from NASA's WMAP satellite [119] 
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bound on the levels of non-gaussinity tnl and (^nl in the trispectrum T^. The amount of 
statistical anisotropy (7^ in the spectrum is in the range g(; ~ 0.290 ± 0.093 [79] . 

Regarding the statistical descriptors, non-gaussianity in the primordial curvature pertur- 
bation ( is one of the subjects of more interest in modern cosmology, because the non- 
gaussianity parameters /nl and tnl together with the spectrum amplitude and spectral 
index allow us to discriminate between the different models proposed for the origin of the 
large-scale structure (see for example Refs. [3],[5l[6]). The most studied and popular models 
are those called the slow-roll models with canonical kinetic terms, because of their simplic- 
ity and because they easily satisfy the spectral index requierements from observation. 
However, the usual predictions of these models is that the levels of non-gaussianity in the 
primordial curvature perturbation are expected to be unobservable [221 1148^ 1191^ 1211^ 1226] . 
However, as we will show in chapters [3] and HI there are some aditional issues that have not 
been taken into account in the current literature. We study these issues to show that it is 
possible to generate sizeable and observable levels of non-gaussianity in a subclass of small- 
field slow-roll inflationary models with canonical kinetic terms; our main conclussion is that 
if non-gaussianity is detected, the aforementioned models could have strong possibilities to 
be the ones responsibles for the formation of the large-scale structure. 

According to the usual assumption, one or more of these scalar field perturbations are 
responsible for the curvature perturbation. In that case, the n-point correlators of ( are 
translationally and rotationally invariants. However, violations of such invariances entail 
modifications of the usual definitions for the spectral functions in terms of the statistical 
descriptors [T1[T2|,|1T]. These violations may be consequences either of the presence of vector 
field perturbations [T2l [M Ii9l 1501 1^1 1521 [5^ 17^ 171 1751 im M 1951 TO 
1224] . spinor field perturbations [301 1194] . or p-form perturbations [TQl [TH lllll I115[ 1116] . 
contributing significantly to (, of anisotropic expansion [TTl [301 El EH EH I102[ I115[ 11651 
11661 1218] or of an inhomogeneous background [121 SH Violation of the statistical 

isotropy (i.e. violation of the rotational invariance in the n-point correlators of () seems 
to be present in the data [HI [8OI [881 1178] and, although its statistical significance is still 
low, the continuous presence of anomalies in every CMB data analysis (see for instance 
Refs. [3l[58l[62l[63l[85l[86l[9ll[92l [1221 [ISSl MM 1205] ) suggests the evidence might be 
decisive in the forthcoming years. The presence of vector fields in the inflationary dynamics 
is not only important to be responsible of violations of the statistical isotropy, they also 
may generate sizeable levels of non-gaussianity described by /nl and tnl; particularly we 
will show in Chapter [5l that including vector fields allows us to get consistency relations 
between the statistical descriptors, more precisely between the non-gaussianity levels /nl 
and Tnl and the amount of statistical anisotropy (7^. 

Because of the progressive improvement in the accuracy of the satellite measurements de- 
scribed above, it is pertinent to study the statistical descriptors of the primordial curvature 
perturbation (, for cosmological models of the origin of the large-scale structure in the 
Universe. It is very important because they could be a crucial tool to discriminate between 
some of most usual cosmological models [SI [S]- 
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The layout of the thesis is the following: the Chapter |2] is devoted to study the statistical 
descriptors for a probability distribution function and its relation with the observational 
parameters, i.e., the spectrum amplitude, the spectral index the levels of non-gaussianity 
/nl and tnl in the bispectrum i?^ and trispectrum respectively and the level of statistical 
anisotropy in the power spectrum, (7^. In this chapter, we also review some generalities of 
the formalism, it has become the standard technique to calculate C and its statistical 
descriptors. In Chapters [3] and H] we show that it is possible to attain very high, including 
observable, values for the levels of non-gaussianity Jml and tnl, in a subclass of small- 
field slow-roll models of inflation with canonical kinetic terms. Comparison with current 
observationally bounds is made. Chapter \5\ is devoted to study the statistical descriptors 
of the primordial curvature perturbation ( when scalar and vector fields perturbations 
are present in the inflationary dynamicc. The levels of non-gaussianity /nl and tnl are 
calculated and related to the level of statisitcal anisotropy in the power spectrum, (7^. We 
show that the levels of non-gaussianity may be very high, in some cases exceeding the 
current observationally limit. Finally we conclude in Chapter [61 



Chapter 



6N FORMALISM AND STATISTICAL 
DESCRIPTORS FOR ( 



I SECTION 2.1 

Introduction 



The primordial curvature perturbation (, as well as the contrast in the temperature of 
the cosmic microwave background radiation 6T/T and the gravitational potential are 
examples of cosmological functions of space and time being described by probability dis- 
tribution functions. In particular, the probability distribution function f{() for ( has well 
defined statistical descriptors which depend directly upon the particular inflationary model 
and that are suitable for comparison with present observational data. Such a comparison 
allows us either to reject or to keep particular inflationary models as those which better 
represent nature's behaviour. In this chapter we give a complete and general description 
of the 5N formalism including both scalar and vector field^l. This formalism is a powerful 
tool and it is commonly used to calculate the primordial curvature perturbation (. We 
also present a cosmologically motivated description of the statistical descriptors for prob- 
ability distribution functions, focusing mainly on /(C)- Finally, in Section l?!4] we give the 
observational constraints for (. 



^ Vector fields will be responsible of violations of the statistical isotropy, as we will see in Chapter [31 
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— I section1!9 

The 6N formalism 



The 5N formalism [521 UMl EM MH [M EDS] provides a powerful method for calculating 
C and all its statistical descriptors at any desired order in cosmological perturbation theory. 
The SN formalism for scalar field perturbations was given at the linear level in Refs. |18H 
1202] and at the non-linear level, which generates non- gaussianity, it was described in Refs. 
|139[ 1142] . In a recent paper |52] the 5N formalism was extended to include vector as well 
as scalar fields. In this section, we give a brief review of the formalism without assuming 
statistical isotropy and define the primordial curvature perturbation (. 

In the cosmological standard model |155] the observable Universe is homogeneous and 
isotropic, being described by the unperturbed Friedmann-Robertson- Walker metric whose 
line element, for a spatially flat universe, looks as follows: 

ds'^ = -dt^ + a\t)6ijdx'dx^ , (2.1) 

where a{t) is the global expansion parameter, t is the cosmic time, and x represents the 
position in cartesian spatial coordinates. The homogeneity and isotropy conditions describe 
very well the Universe at large scales, but departures from the unperturbed background are 
observationally evident at smaller scales. 

One way to parametrize the departures from the homogeneous and isotropic background is 
to include perturbations in the metric, for which we have to define a slicing and a threading. 
The slicing will be defined so that the energy density in fixed-t slices of spacetime is uniform. 
The threading will correspond to comoving fixed-x world lines. With generic coordinates 
the perturbed metric of the perturbed universe may be defined as 

ds'^ = g^ydx^dx" . (2.2) 

To define the cosmological perturbations, one chooses a coordinate system in the perturbed 
universe, and then compares that universe with an unperturbed one. The unperturbed 
universe is taken to be homogeneous, and is usually taken to be isotropic as well. 

p. 2.1 The curvature perturbation 

To define the curvature perturbation, we smooth|§the metric tensor and the energy- momentum 
tensor on a comoving scale R and one considers the super- horizon regime aR > H^^ where 



^Smoothing a function ^(x) means that g at each location is replaced by its average within a sphere 
of coordinate radius R around that position. The averaging may be done with a smooth window function 



The SN formalism 



10 



H = a/ais the Hubble parameter and a{t) is the scale factor normalised to 1 at presenlH. 
The energy density p and pressure P are smoothed on the same scale. On the reason- 
able assumption that the smoothing scale is the biggest relevant scale, spatial gradients 
of the smoothed metric and energy-momentum tensors will be negligible. As a result, the 
evolution of these quantities at each comoving location will be that of some homogeneous 
'separate universe'. In contrast with earlier works on the separation universe assumption, 
we will in this section allow the possibility that the separate universes are anisotropic even 
though homogoneous. 

We consider the slicing of spacetime with uniform energy density, and the threading which 
moves with the expansion (comoving threading). By virtue of the separate universe as- 
sumption, the threading will be orthogonal to the slicing. The spatial metric can then be 
written as 

^.,(x,t) = a^(t)e^«-*)(Je^'^(-*))^^., (2.3) 

where I is the unit matrix, and the matrix h is traceless so that le^^ has unit determinant. 
The smoothing scale is chosen to be somewhat shorter than the scales of interest, so that 
the Fourier components of C, on those scales is unaffected by the smoothing. The time 
dependence of the locally defined scale factor a(x, t) = a{t) exp(C) defines the rate at which 
an infinitesimal comoving volume V expands: V/V = 3a(x, t)/a(x, t). 

We split In a and hij into an unperturbed part plus a perturbation: 

In a(x,r) = In a(r) + C(x, r) , (2.4) 
hij{yi,T) = hij{T) ^ 5hij{yi,T) . (2.5) 

The unperturbed parts can be defined as spatial averages within the observable Universe, 
but any definition will do as long as it makes the perturbations small within the observable 
Universe. If they are small enough, ^ and 5hij can be treated as first-order perturbations. 
That is expected to be the case, with the proviso that a second-order treatment of C, will 
be necessary to handle its non-gaussianity if that is present at a level corresponding to 
/nl ^ 1 (with the gaussian and non-gaussian components correlated) |143] . 

Under the reasonable assumption that the Hubble scale is the biggest relevant distance 
scale, the energy continuity equation diVp) = —PdV at each location is the same as in a 
homogeneous universe; as far as the evolution of p is concerned, we are dealing with a 
family of separate homogeneous universes. With the additional assumption that the initial 
condition is set by scalar fields during inflation, the smoothed hij{t) is time-independent 
after smoothing and then the separate universes are homogeneous as well as isotropic. 



such as a gaussian. The smoothed function is supposed to have no significant Fourier components with 
coordinate wavenumbers k 3> R~^, which means that its gradient at a typical location is at most of order 
1/R. A function g with that property is said to be 'smooth on the scale R'. 
^The value of the scale factor at present epoch is usually denoted by ao- 
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Since we are working on slices of uniform p, the energy continuity equation can be written 

p{t) = -3 [Hit) + C(x, t)] [pit) + P(x, t)] . (2.6) 

One write 

P(x,t)=P(t) + 5P(x,t), (2.7) 

so that SP is the pressure perturbation on the uniform density shces, and choose P(t) so 
that the unperturbed quantities satisfy the unperturbed equation p = —3Hip + P). Then 

• H5P , , 

This gives C if we know pit) and P(x, t). It makes C(x) time-independent0 during any 
era when P is a unique function of p |139l 1219] (hence uniform on slices of uniform p). 
The pressure perturbation is said to be adiabatic in this case, otherwise it is said to be 
non-adiabatic. 



The key assumption in the above discussion is that in the superhorizon regime certain 
smoothed quantities (in this case p and P) evolve at each location as they would in an 
unperturbed universe. In other words, the evolution of the perturbed universe is that of a 
family of unperturbed universes. This is the separate universe assumption, that is useful 
also in other situations [140j . 

The primordial curvature perturbation ^ is directly probed by observation on 'cosmological 
scales' corresponding to roughly e'^^H^^ < < -ffo~1fl- These scale begin to enter the 
horizon when when T ~ 1 MeV. The Universe at that stage is radiation dominated to very 
high accuracy [m [991 1107[ 1108] . implying P = p/3 and a constant curvature perturbation 
which we denote simply by C(x), and it is the one constrained by observation as we will see 
in section 12.41 When cosmological scales are the only ones of interest, one should choose 
the smoothing scale as P ~ c'^^Hq^. Unless stated otherwise, we make this choice. 

Within a given scenario, ( will exist also on smaller inverse wavenumbers, down to some 
'coherence length' which might be as low as ~ c'^^Hq^ (the scale leaving the horizon 
at the end of inflation). If one is interested in such scales, the smoothing scale R should be 
chosen to be (somewhat less than) the coherence length. 



f .2.2 The SN formula 

Keeping the comoving threading, we can write the analogue of Eq. (12. 3p for a generic slicing: 

^.,(x,r)^a2(x,r)(/e2^(-'-))^ , (2.9) 



'^Absorving ( into the unperturbed scale factor. 
^Hq is the Hubble parameter today. 
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with again le^^ having unit determinant so that the rate of volume expansion is V/V = 
3a(x, t). Starting with an initial 'flat' slicing such that the locally-defined scale factor is 
homogeneous, and ending with a slicing of uniform density, we then have 

C(x, t) = 5iV(x, t) = iV(x, t) - iVo(t), (2.10) 

where the number of e-folds of expansion is defined in terms of the volume expansion by 
the usual expression N = V/3V. The choice of the initial epoch has no effect on 6N, 
because the expansion going from one fiat slice to another is uniform. We will choose the 
initial epoch to be a few Hubble times after the smoothing scale leaves the horizon during 
inflation. According to the usual assumption, the evolution of the local expansion rate is 
determined by the initial values of one or more of the perturbed scalar fields Then we 
can write 

M^) = (f>i + SM^), (2.11) 

C(x,t) = 5iV(0i(x),02(x),...,t) = iV/(t)50/(x) + 

+ ^Njj{t)6M^)SM^) + -^NjjK{t)6M^)SM^)SM^) • • • , (2.12) 

where Nj = dN/dcpi, etc., and the partial derivatives are evaluated with the fields at their 
unperturbed values denoted simply by </>/. The field perturbations 6 (pi in Eq. (12.12^ are 
defined on the 'flat' slicing such that a(x, t) is uniform. 

The unperturbed fleld values are deflned as the spatial averages, over a comoving box 
within which the perturbations are deflned. The box size aL should satisfy LHq ^ 1 
so that the observable Universe should flt comfortably inside it |138] . Observations are 
available within the observable universe and, except for the low multipoles of the CMB, 
all observations probe scales k ^ Hq. To handle them, one should choose the box size as 
L = Hq^ |121] . a smaller choice would throw away some of the data while a bigger choice 
would make the spatial averages unobservable. Low multipoles i of the CMB anisotropy 
explore scales of order Hq^/E not very much smaller than H^^. To handle them one has to 
take L bigger than Hq^. For most purposes, one should use a box, such that In(Lifo) is just 
a few (ie. not exponentially large) [1141 1138^ I14U] . When comparing the loop contribution 
with observation one should normally set L = Hq^, except for the low CMB multipoles 
where one should choose L ^ Hq^ with ln(/cL) ~ 1. With the choice L = Hq^, ln(/cL) ~ 5 
for the scales explored by the CMB multipoles with i ~ 100, while In(fcL) ~ 10 for the 
scales explored by galaxy surveys. Since we are interested in giving orders of magnitude 
and simple mathematical expressions, in the current thesis we will set ln(A;L) ~ 10 without 
loss of generality. 

The spatial averages of the scalar flelds, that determine Nj, etc., and hence ( cannot in 
general be calculated. Instead they are parameters, that have to be specifled along with 
the relevant parameters of the action before the correlators of ( can be calculated. The 



^As we will see in the next chapters, the ln(A:L) dependence appear when we consider higher order 
corrections to C, explicitly in the higher order correlators which describe its statistical properties. 
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only exception is when ( is determined by the perturbation of the inflaton in single-field 
infiation. Then, the unperturbed field value when cosmological scales leave the horizon can 
be calculated, knowing the number of e-folds to the end of infiation which is determined 
by the evolution of the scale factor after infiation. Although the unperturbed field values 
cannot be calculated, their mean square for a random location of the minimal box (ie. of 
the observable Universe) can sometimes be calculated using the stochastic formalism 



On the other hand, if we suppose that one or more perturbed vector fields also affect the 
evolution of the local expansion rate, the curvature perturbation, in the simplest case where 
( is generated by one scalar field and one vector field and assuming that the anisotropy 
in the expansion of the Universe is negligible, can be calculated up to quadratic terms by 
means of the following truncated expansion 



C(x) = 5iV(0(x),A,(x),t) 

+ N'JA, + lN^^{6(j)f + NiJ<f)6A, + ^-N^i^dAdA, , (2.13) 



where 



del) ' ^ " dAi ' '^'^ ~ ' ''^^ - dAidA, ' " dAidi 



(j) being the scalar field and A the vector field, with i denoting the spatial indices running 
from 1 to 3. As with the scalar fields, the unperturbed vector field values are defined as 
averages within the chosen box. 

In these formulas there is no need to define the basis (triad) for the components Ai. Also, 
we need not assume that A^ comes from a 4-vector field, still less from a gauge field. 



P?5TT The growth of ( 

As noted earlier, ( is constant during any era when pressure P is a unique function of 
energy density p. In the simplest scenario, the field whose perturbation generates ( is the 
infiaton field in a single-field model. Then the local value of is supposed to determine 
the subsequent evolution of both pressure and energy density, making ( constant from the 
beginning. 

Alternatives to the simplest scenario generate all or part of C at successively later eras. 
Such generation is possible during any era, unless there is sufficiently complete matter 
domination (P = 0) or radiation domination (p = P/3). Possibilities in chronological order 
include generation during 

(i) multi-field infiation [202] . 

(ii) at the end of infiation |135j . 

(iii) during preheating [114] . 
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(iv) at reheating, and 

(v) at a second reheating through the curvaton mechanism [1321 11451 1146[ 1153] . 

A vector field cannot replace the scalar field in the simplest scenario, because unperturbed 
infiation with a single unperturbed vector field will be very anisotropic and so will be the 
resulting curvature perturbation. Even with isotropic infiation, we will see in Chapter [5] that 
a single vector field perturbation cannot be responsible for the entire curvature perturbation 
(at least in the scenarios that were discussed in the Ref. [52]) because its contribution is 
highly anisotropic. It could instead be responsible for part of the curvature perturbation, 
through any of the mechanisms listed above. 

I SECTION 2.3j 

Statistical descriptors for a probability distribution function 



A probability distribution function f{() for any function of space and time C(X)^) may be 
understood as the univocal correspondence between the possible values that ( may take 
throughout the space and the normalised frecuency of appearences of such values for a given 
time. Any continous function of ( might represent a probability distribution function as 
long as /(C) > and /(C)'^C = 1- However, for a particular probability distribution 
function, how many independent parameters do we need to completely characterize it in a 
unique way? And despite the possible infinite number of parameters required to do this, 
what is the information encoded in those parameters? 

The answers to these questions rely on the moments m(^{n) of the distribution. 

For a given probability distribution function /(C), there are an infinite number of moments 
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that work as statistical descriptors of ({'x,t): 

the mean value : m^(l) = {0 = JcfiC)dC, (2.15) 

the variance : m^(2) = j {C - iOYfiOdC , (2.16) 



lic -iorm dc 

K(2)]' 



the skewness : m^(3) = , (2-17) 



thekurtosis: m^(4) = (C))'/(C)c?C ^ ^2.18) 



and so on. 

What can we say about C{'^,t) from the knowledge of the moments of the distribution? If, 
for instance, all the odd moments with n > 3 (skewness, ... etc) are zero, we can say that 
the probability distribution function f{() is even around the mean value. If in addition 
all the even moments with n > 4 (kurtosis, ... etc) are expressed only as products of the 
variance, we can say that the distribution function is gaussian. Indeed, as is well known, 
the only quantities required to reproduce a gaussian function are the mean value and the 
variance: 

f . (C) - ^ ^-(C-mc(l))V2mc(2) (2 ]Q) 

A/27rm^(2) 

Departures from the exact gaussianity come either from non-vanishing odd moments with 
n > 3, in which case the probability distribution function is non-symmetric around the mean 
value, or from higher n > 4 even moments different to products of the variance, in which 
case the probability distribution function continues to be symmetric around the mean value 
although its "peakedness"[l| is bigger than that for a gaussian function, or from both of them. 
A non-gaussian probability distribution function requires then more moments, other than 
the mean value and the variance, to be completely reconstructed. Such a reconstruction 
process is described for instance in Ref. |183j . 

Working in momentum space is especially useful in cosmology because the modes associated 
with the quantum fluctuations of scalar flelds during inflation become classical once they 
leave the horizon [3 [Ml EHl El M [1331 [I3S1 [HDl [IM [M [ISD]. The same applies for 
the primordial curvature perturbation ( which, in addition, is a conserved quantity while 
staying outside the horizon if the adiabatic condition is satisfied |139] . As regards the 
moments of the probability distribution function, they have a direct connection with the 
correlation functions for the Fourier modes Ck = / d^k({x.)e~'^^'^ defined in fiat space. 



''Higher even standarized moments different to products of tlie variance mean more of tlie variance is 
due to infrequent extreme deviations, as opposed to frequent modestly- sized deviations. 
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As the ra-point correlators of Ck are generically defined in terms of spectral functions of the 
wavevectors involved: 

two — point correlation — > 

(CkiCka) = 

three — point correlation — )■ 

(CkiCkaCks) = 

four — point correlation — )■ 

(CkiCk2Ck3Ck4) = 



spectrum : 

(2^)353(ki+k2)Pc(k) 
bispectrum 



(2.20) 



i2nf6\k^ + ka + k3)Pc(ki, ka, kg) , (2.21) 
trispectrum : 

(27r)353(ki + k2 + kg + k4)r^(ki, k2, kg, k4) , (2.22) 



and so on, 

the moments of the distribution are then written as momentum integrals of the spectral 
functions for the modes 



the variance : 7ti^(2) 
the skewness : m(^[3) 

the kurtosis : m^(4) 



n(k) 



(27r)3 

J d3fc,g^fe,d3fc3 y^(k^^k2,k3,k4) 



[/0n(k)]' 



(2.23) 
(2.24) 

(2.25) 



and so on. 

Non-gaussianity in ( is, therefore, associated with non-vanishing higher order spectral func- 
tions, starting from the bispectrum B^. 



p. 3.1 Statsitical descriptors for primordial curvature 
perturbation ( 

Theoretical cosmologists work with (. However, astronomers work with observable quanti- 
ties such as the contrast in the temperature of the cosmic microwave background radiation 
6T/T. The connection between the theoretical cosmologist quantity ( and the astronomer 
quantity 6T/T is given by the Sachs- Wolfe effect |176j which, at first-order and for super- 
horizon scales, looks as follows: 




(2.26) 
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Thus, although it is essential to study the Sachs- Wolfe relation at higher orders, which is far 
more complicated than Eq. fl2.26p . theoretical cosmologists may study the statistical prop- 
erties of the observed 6T/T through the spectral functions associated with the curvature 
perturbation 

mean value of 6T/T = — )■ mean value of C = , (2.27) 
variance : msT/T{'^) -> spectrum : (k) , (2.28) 
skewness : msT/ri'^) bispectrum : i?^(ki, k2, ks) , (2.29) 

kurtosis : m5r/r(4) — )■ trispectrum : T^(ki, k2, ka, k4) , (2.30) 



and so on. 

Now, we will parametrize the spectral functions of ( in terms of quantities which are the 
ones for which observational bounds are given. Because of the direct connection between 
these quantities and the moments of the probability distribution function /(C), we may also 
call these quantities as the statistical descriptors for /(C)- 

The bispectrum and trispectrum are parametrized in terms of products of the spec- 
trum P^, and the quantities /nl and tnl and g^^i^ respectivelj|§ [HTl HUI [1^5] : 

P,<ki, k2, ks) = ^/jVL(kl, k2, ks) [Pc(ki)Pc(k2) + c. p.] , (2.31) 
T^(ki, k2, ks, k4) = ^rjVL(kl, k2, ks, k4) [Pc(ki)P^(k2)Pc(ki + k4) + c. p.] + 

54 

+ -^77VL(ki,k2,k3,k4) [Pc(ki)Pc(k2)Pc(ks) +C. p.] , 

(2.32) 

where c. p. means cyclic permutations. Higher order spectral functions would be parametrized 
in an analogous way. 



By virtue of the reality condition C(~k) = C*(k), an equivalent definition of the spectrum 
is 



(C(ki)C*(k2)) = {2n)'6'{k^ - k2)Pc(k) . (2.33) 



^If we assume that the statistical inhomogeneity is present, i.e translational invariance of the n-point 
correlators of ^ is broken, it is necesary introduce modifications of the usual definitions of the statistical 
descriptors of the primordial curvature perturbation ^. For example, if there is statistical inhomogeneity 
then (CkiCka) is not proportional to 6{ki + k2). In this thesis we will assume statistical homogeneity, so 
that the statistical descriptors given in the last subsection can be correctly applied. 

^There is actually a sign difference between the /nl defined here and that defined in Ref. |148j . The 
origin of the sign difference lies in the way the observed Jnl is defined |I20j . through the Bardeen's 
curvature perturbation: = $f + ff{L{^fY with = (3/5)C, and the way Jnl is defined in Ref. 
[148], through the gauge invariant Newtonian potential: = $^ + JnU^l)^ with $^ = -(3/5)C. 
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Setting ki = k2 the left hand side is (|C(k)|^). It follows that the the spectrum is positive 
and nonzero. Even if -P(^(k) is anisotropic, the reality condition requires -P((k) = P(^(— k). 
The spectrum will therefore be of the form [1] 



k and g(; is the level of statistical anisotropy. The homogeneity and isotropy requirements at 
large scales imply that the spectrum and bispectrum are functions of the wavenumbers 
only. For the trispectrum and the other higher order spectral functions, the momentum 
dependence also involves the direction of the wavevectors. 

On the other hand if the n-point correlators are also invariant under rotations (statistical 
isotropy) the spectral functions Pf(k) = (27r^//c^)P^(k) and i?^(ki, k2, ka) depend only on 
the magnitude of the wavevectors [1] . In this case the spectrum is parametrized in terms 

1/2 

of an amplitude and a spectral index which measures the deviation from an exactly 
scale- invariant spectrum |14Uj : 



Given the present observational status, n^, Jnl, tnl, Qnl and are the statistical descrip- 

1/2 

tors that discriminate among models for the origin of the large-scale structure once 
has been fixed to the observed value. Since non-vanishing higher order spectral functions 
such as Bi^ and imply non-gaussianity in the primordial curvature perturbation (, the 
statistical descriptors /nl, tnl and (^nl are usually called the levels of non-gaussianity. 

^ SECTION 2.4 , 

Observational constraints 



Direct observation, coming from the anisotropy of the CMB and the inhomogeneity of 
the galaxy distribution, gives information on what are called cosmological scales |119j . 
These correspond to a range A In A; ~ 10 or so downwards from the scale k"^ ~ Hq^ that 
corresponds to the size of the observable Universe. 





(2.35) 




Spectrum and non-gaussianity 



Observational results concerning the spectrum V( are generally obtained with the assump- 
tion of statistical isotropy ((7^ = 0), but they would not be greatly affected by the inclusion 
of anisotropy at the 10% level. 
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NASA's COBE-satellite provided us with a reliable value for the spectral amplitude 
[35]: ^ (4.957 ± 0.094) x 10"^ which is usually called the COBE normalisation. As 

regards the spectral index, the latest data release and analysis from the WMAP satellite 
shows that = 0.960 ± 0.014 |119] which rejects exact scale invariance at more than 2a. 
Such a result has been extensively used to constrain inflation model building [5], and al- 
though several classes of inflationary models have been ruled out through the spectral index, 
lots of models are still allowed; that is why it is so important an appropiate knowledge of 
the statistical descriptors Jnl and t^l. Present observations show that the primordial cur- 
vature perturbation ( is almost, but not completely, gaussian. The level of non-gaussianity 
/jvL in the bispectrum 5^, after five years of data from NASA's WMAP satellite, is in the 
range —9 < Jnl < HI at 2a |119j . There is at present no observational bound on the 
level of non-gaussianity tnl in the trispectrum although it was predicted that COBE 
should either detect it or impose the lower bound Ittv^,! < 10* [3T | I163j . It is expected that 
future WMAP data releases will either detect non-gaussianity or reduce the bounds on f^i 
and tnl at the 2a level to I/^vlI < 40 [I2D] and It^vlI < 2 x 10^ [112] respectively. The 
ESA's PLANCK satellite [23 1206] , launched in 2009, promises to reduce the bounds to 
I/atlI < 10 [120j and Itat^^I < 560 [112] at the 2a level if non- gaussianity is not detected. 
In addition, by studying the 21-cm emission spectral line in the cosmic neutral Hydrogen 
prior to the era of reionization, it is also possible to know about the levels of non-gaussianity 
/tvl and tnl, the 21-cm background anisotropies capture information about the primordial 
non-gaussianity better than any high resolution map of cosmic microwave background ra- 
diation: an experiment like this could reduce the bounds on the non-gaussianity levels to 
IInl] ^ 0.2 [m Us] and \tnl\ ^ 20 [15] at the 2a confidence. Finally, it is worth stating 
that there have been recent claims about the detection of non-gaussianity in the bispectrum 
of ( from the WMAP 3-year data [1011 1223j . Such claims, which report a rejection of 
Inl = at more that 2a (26.9 < fj^L < 146.7), are based on the estimation of the bispec- 
trum while using some specific foreground masks. The WMAP 5-year analysis [119j shows 
a similar behaviour when using those masks, but reduces the significance of the results 
when other more conservative masks are included allowing again the possibility of exact 
gaussianity. 




Statistical anisotropy and statistical inhomogeneity 



Taking all the uncertainties into account, observation is consistent with statistical anisotropy 
and statistical inhomogeneity allowing either of these things at around the 10% level. In- 
deed, some recent papers [HI [79], [HQ], [HHl I178j claim for the presence of statistical anisotropy 
in the five-year data from the NASA's WMAP satellite [159j . In this section we briefly re- 
view what is known. 

Assuming statistical homogeneity of the curvature perturbation, a recent study [79] of the 
CMB temperature perturbation finds weak evidence for statistical anisotropy (see Fig. 12.11 
for an example of who is statitiscal anisotropy). They keep only the leading (quadrupolar) 
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(a) (b) (c) 

Figure 2.1: Statistical anisotropy. (a). This simulation corresponds to the statistical 
isotropic case = 0. (b). In contrast, this simulation corresponds to the statistical 
anisotropic case with g^^ = 1, d pointing along the horizontal direction and setting to zero 
the isotropic part. (c). Same as in (b) but with d pointing along the vertical direction. 
(Courtesy of Mindaugas Karciauskas) . 



term of Eq. (ICTj) : 

Pc(k)=Pf(A:)(l + (7c(d-k)2) , (2.36) 

and find g(; ~ 0.290 ±0.031 which rules out statistical isotropy at more than 9a. Neverthe- 
less, the preferred direction lies near the plane of the solar system, which makes the authors 
of Ref. [79] believe that this effect could be due to an unresolved systematic error (among 
other possible systematic errors which have not been demonstrated either to be the source 
of this statistical anisotropy nor to be completely uncorrelated [7^). 

Even if the result found in Ref. [79] turns out to be due to a systematic error, some fore- 
casted constraints on show that the statistical anisotropy subject is worth studying |167] : 
\g(;\ < 0.1 for the NASA's WMAP satellite [159] if there is no detection, and l^f^l < 0.02 for 
the ESA's PLANCK satellite [52] if there is no detection. There is at present no bound on 
statistical anisotropy of the 3-point or higher correlators. 

In some different studies, the mean-square CMB perturbation in opposite hemispheres has 
been measured, to see if there is any difference between hemispheres. Quite recent works 
[62| [63| |86] find a difference of order ten percent, for a certain choice of the hemispheres, 
with statistical significance at the 99% level. Given the difficulty of handling systematic 
uncertainties it would be premature to regard the evidence for this hemispherical anisotropy 
as completely overwhelming. Nevertheless, what would hemispherical anisotropy imply for 
the curvature perturbation? Focussing on a small patch of sky, the statistical anisotropy of 
the curvature perturbation implies that the mean-square temperature perturbation within 
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a given small patch will in general depend on the direction of that patch. This is because 
the mean square within such a patch depends upon the mean square of the curvature 
perturbation in a small planar region of space perpendicular to the line of sight located at 
last scattering. But the mean-square temperature will be the same in patches at opposite 
directions in the sky, because they explore the curvature perturbation C(k) in the same 
k-plane and the spectrum V(^{k.) is invariant under the change k — )■ — k. It follows that 
statistical anisotropy of the curvature perturbation cannot by itself generate a hemispherical 
anisotropy. We may then conclude that hemispherical anisotropy of the CMB temperature 
requires statistical inhomogeneity of the curvature perturbation. Then (C(ki)C(k2)) is not 
proportional to 5'^(ki + k2) [H]. 



I SECTION 2.5^ 



Conclusions 



In this chapter we have presented the statistical description of primordial curvature per- 
turbation (. This framework allows us to define statistical descriptors, which provides a 
bridge between theory and observation. The relevant parameters that parametrize these 
statistical descritors are: the spectrum amplitude Vq^^'^ , the spectral index n^, the levels of 
non-gaussianity /nl, tnl and (^nl and the level of statistical anisotropy g(^. Some of these 
parameters have an upper bound from observation, so it is very important to study the- 
oretical models that successfully reproduce these observations. To study these theoretical 
aspects, we have presented a poweful tool to calculate the primordial cuvature pertubation 
C, and all its statistical descriptors; such a tool is usuUy called the 5N formalism. In the 
next two chapters we will see how to work out this formalism. 



Chapter 



ON THE ISSUE OF THE ( SERIES 
CONVERGENCE AND LOOP CORRECTIONS IN 
THE GENERATION OF OBSERVABLE 
PRIMORDIAL NON- GAUSSIANITY IN 
SLOW-ROLL INFLATION: THE BISPECTRUM 



SECTION 3.1 

Introduction 



Since COBE [158] discovered and mapped the anisotropics in the temperature of the cosmic 
microwave background radiation |199j . many balloon and satellite experiments have refined 
the measurements of such anisotropics, reaching up to now an amazing combined precision. 
The COBE sequel has continued with the WMAP satellite |159] which has been able to 
measure the temperature angular power spectrum up to the third peak with unprecedent 
precision ^U\, and increase the level of sensitivity to primordial non-gaussianity in the 
bispectrum by two orders of magnitude compared to COBE [118(1119] . The next-to-WMAP 
satellite, PLANCK [59], which was launched in may of 2009, is expected to precisely measure 
the temperature angular power spectrum up to the eighth peak |206j . and improve the level 
of sensitivity to primordial non-gaussianity in the bispectrum by one order of magnitude 
compared to WMAP [HP] . 
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Because of the progressive improvement in the accuracy of the sateUite measurements de- 
scribed above, it is pertinent to study cosmological inflationary models that generate signif- 
icant (and observable) levels of non-gaussianity. An interesting way to address the problem 
involves the SN formalism [521 1139[ I142[ I18H I182[ 1202] . which can be employed to give the 
levels of non-gaussianity /tvl |142] and TAr^ EI] in the bispectrum and trispectrum 
of the primordial curvature perturbation ( respectively. Such non-gaussianity levels are 
given, for slow-roll inflationary models, in terms of the local evolution of the universe un- 
der consideration, as well as of the n-point correlators, evaluated a few Hubble times after 
horizon exit, of the perturbations S(j)i in the scalar fields that determine the dynamics of 
such a universe during inflation. 

In the 5 A'^ formalism for slow-roll inflationary models, the primordial curvature perturbation 
C(x, t) is written as a Taylor series in the scalar field perturbations (50i(x, t^,) evaluated a 
few Hublle times after horizon exit0, 

IJ IJ 

■ UK 

- ^^iV/jx(t)(<50/(x,t^)50j(x,t,)50i^(x,t^)) + 

■ UK 

+... , (3.1) 

It is in this way that the correlation functions of ( (for instance, (CkiCk2Ck3)) can be ob- 
tained in terms of series, as often happens in Quantum Field Theory where the probability 
amplitude is a series whose possible truncation at any desired order is determined by the 
coupling constants of the theory. A highly relevant question is that of whether the series for 
6N converges in cosmological perturbation theory and whether it is possible in addition to 
find some quantities that determine the possible truncation of the series, which in this sense 
would be analogous to the coupling constants in Quantum Field Theory. In general such 
quantities will depend on the specific inflationary model; the series then cannot be simply 
truncated at some order until one is sure that it does indeed converge, and besides, one has 
to be careful not to forget any term that may be leading in the series even if it is of higher or- 
der in the coupling constant. This issue has not been investigated in the present hterature, 
and generally the series has been truncated to second- or third-order neglecting in addition 
terms that could be the leading ones [3 El [221 [3ll SQl [US [Ml MSi [HIl [2251 [2261 1228] . 

The most studied and popular inflationary models nowadays are those of the slow-roll vari- 
ety with canonical kinetic terms |137[ I140[ 1141] , because of their simplicity and because they 



^This equation is similar to one in Eq. (|2.12p . the difference is tliat fiere we are redefined it so tfiat 
(C(t,x)=0). 
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easily satisfy the spectral index requirements for the generation of large-scale structures. 
One of the usual predictions from inflation and the theory of cosmological perturbations is 
that the levels of non-gaussianity in the primordial perturbations are expected to be unob- 
servably small when considering this class of models [221 EM CHHl HHHl HHQl IMl [ml 
I226[ I228j^ . This fact leads us to analyze the cosmological perturbations in the framework 
of first-order cosmological perturbation theory. Non-gaussian characteristics are then sup- 
pressed since the non-linearities in the infiaton potential and in the metric perturbations 
are not taken into account. The non-gaussian characteristics are actually present and they 
are made explicit if second-order |143j or higher-order corrections are considered. 

The whole literature that encompasses the slow-roll inflationary models with canonical 
kinetic terms reports that the non- gaussianity level f^L is expected to be very small, 
being of the order of the slow- roll parameters and rji, (e^, \rii\ -C 1) [22 | I148[ [19H I21H 1226] . 
These works have not taken into account either the convergence of the series for ( nor the 
possibility that loop corrections dominate over the tree level ones in the n-point correlators. 
Our main result in this chapter is the recognition of the possible convergence of the ( series, 
and the existence of some "coupling constants" that determine the possible truncation of 
the ( series at any desired order. When this situation is encountered in a subclass of small- 
field slow-roll inflationary models with canonical kinetic terms, the one-loop corrections 
may dominate the series when calculating either the spectrum P<^, or the bispectrum 5^. 
This in turn may generate sizeable and observable levels of non-gaussianity in total contrast 
with the general claims found in the present literature. 

The layout of the chapter is the following: Section 13.21 is devoted to the issue of the ( series 
convergence and loop corrections in the framework of the 6N formalism. The presentation 
of the current knowledge about primordial non-gaussianity in slow-roll inflationary models 
is given in Section [3l3l A particular subclass of small-field slow- roll inflationary models is 
the subject of Section [23] as it is this subclass of models that generate significant levels of 
non- gaussianity. The available parameter space for this subclass of models is constrained 
in Section 13.5! by taking into account some observational requirements such as the COBE 
normalisation, the scalar spectral tilt, and the minimal amount of inflation. Another re- 
quirement of methodological nature, the possible tree-level or one-loop dominance in 
and/or P^, is considered in this section. The level of non-gaussianity /tvl in the bispec- 
trum is calculated in Section 13.61 for models where ( is generated during inflation; a 
comparison with the current literature is made. Section 13.7! is devoted to central issues 
in the consistency of the approach followed such as satisfying necessary conditions for the 



^One possible exception is the two-field slow-roll model analyzed in Ref. |5] (see also Refs. |28l[29]) where 
observable, of order one, values for fwL are generated for a reduced window parameter associated with the 
initial field values when taking into account only the tree- level terms in both P(^ and B(^. However, such a 
result seems to be incompatible with the general expectation, proved in Ref. |211| . of f^L being of order 
the slow-roll parameters, and in consequence unobservable, for two-field slow-roll models with separable 
potential when considering only the tree- level terms both in and B(^. The origin of the discrepancy 
could be understood by conjecturing that the trajectory in field space, for the models in Refs. |SJ [IHl [H], 
seems to be sharply curved, being quite near a saddle point; such a condition is required, according to Ref. 
Uni, to generate /nl ^ 0{l). 
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convergence of the ( series and working in a perturbative regime. Finally in Section 13.81 
we conclude. As regards the level of non-gaussianity tnl in the trispectrum T^, it will be 
studied in the following Chapter. 



SECTION 3.2 



( series convergence and loop corrections 



In order to calculate ({t, x) from Eq. fl2.12p . we need information about the physical content 
of the Universe at times t and tin- By choosing the initial time tin a few Hubble times after 
the cosmologically relevant scales leave the horizon during inflation tin = t*, and the final 
time t corresponding to a slice of uniform energy density, we recognize that A^, for slow-roll 
inflationary models, is completely parametrized by the values a few Hubble times after 
horizon exit of the scalar fields (pi present during inflation and the energy density at the 
time at which one wishes to calculate (: 

C{t, x) = N{p{t), Mt*, x), 02(t., x), ...) - N{p{t), MQ, MQ, •••) • (3.2) 

The previous expression can be Taylor-expanded around the unperturbed background val- 
ues for the scalar fields (pi and suitably redefined so that {({t,x.)) = 0. Thus, 

C(t,x) = J2Ni{t)6(Pi{^,Q-J2Ni{t){6(Pi{^,Q) + 
I I 

+^J2Nij{t)6(Pj{^,Q6(Pj{^,Q -^J2Nu{t){6(Pi{^,Q6(Pj{^,Q) + 
IJ IJ 

+^ J2 NijK{t)6(Pi{^, U)6(Pj{^, QS(Pk{^, Q 

■ UK 

- ^$^A'wW(50/(x,t*)50j(x,tJ50i^(x,t,)) + 

■ UK 

+... , (3.3) 

where the 6(pi{t^,,x.) are the scalar field perturbations in the flat slice a few Hubble times 
after horizon exit, whose spectrum amplitude is given by [33] 

< = (3.4) 

and the notation for the derivatives is A^,- = -f^. A,-, = aTaj. i and so on. 



The expression in Eq. (13. 3p has been used to calculate the statistical descriptors of C at 
any desired order in cosmological perturbation theory by consistently truncating the series 
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|142j . For instance, by truncating the series at first order, the amphtude of the spectrum 
of C defined in Eqs. fl2:20|) and fl2:35|) is given by [181] 

(3.5) 

which in turn gives the well known formula for the spectral index [181] : 

n,-l = -2e- ^^P yj^/^f ' (3-6) 

where a subindex i in means a derivative with respect to the (pi field, and being e one of 
the slow-roll parameters defined by e = —H/H'^, mp = (87rG)~^ the reduced Planck mass, 
and V the scalar inflationary potential. Analogously, the level of non-gaussianity f^L in the 
bispectrum of ( defined in Eqs. ( I2.2ip and fl2.3ip is obtained by truncating the series 
at second order and assuming that the scalar field perturbations 6(f)i are perfectly gaussian 

m- 

In the last expression the ln(A;L) factor is of order one, L being the infrared cutoff when 
calculating the stochastic properties in a minimal box [271 1138] . 

The truncated series methodology has proved to be powerful and reliable at reproducing 
successfully the level of non-gaussianity f^i in single-field slow-roll models |189] and in the 
curvaton scenario [31]. Nevertheless, for more general models, how reliable is it to truncate 
the series at some order? In the first place, from Eq. fl3.3p it is impossible to know whether 
the series converges until the N derivatives are explicitly calculated and the convergence 
radius is obtained; obviously if the series is not convergent at all, the expansion in Eq. 
(13. 3p is meaningless. Without any proof of the contrary, the current assumption in the 
literature [3 El [221 [SH SQl Iffl MM [US 12111 [2251 [228] has been that the C series 
is convergent. In addition, supposing that the convergence radius is finally known, the 
truncation at any desired order would again be meaningless if some leading terms in the 
series get excluded. Such a situation might easily happen if each x-dependent term in the ( 
series is considered smaller than the previous one, which indeed is the standard assumption 
[3 [61 [221 [an [101 [HD [Ha mm MSI MS\ , but which is not a universal fact. 

When studying the series through a diagrammatic approach [39], in an analogous way 
to that for Quantum Field Theory via Feynman diagrams, the first-order terms in the 
spectral functions are called the tree-level terms. Examples of these tree-level terms are 
those in Eqs. (13. 5 p and (13. 6p . and the first one in Eq. (13. 7p . Higher-order corrections, 
such as that which contributes with the second term in Eq. (13. 7p . are called the loop 
terms because they involve internal momentum integrations. The statistical descriptors of 
( has been so far studied by naively neglecting the loop corrections against the tree-level 
terms [SI [SI [221 [SB [IDl MH [M MM [223 [22S1 [22B]; nevertheless, as might happen 
in Quantum Field Theory, eventually some loop corrections could be bigger than the tree- 
level terms, so it is essential to properly study the possible n-loop dominance in the spectral 
functions. 
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I SECTION 



Non-gaussianity in slow-roll inflation 



The most frequent class of inflationary models found in the literature are those which satisfy 
the so called slow-roll conditions, as these very simple models easily meet the spectral index 
observational requirements discussed in Subsection 12.41 for the generation of large-scale 
structures. 

The slow-roll conditions for single-field inflationary models with canonical kinetic terms 
read 

0' « V{<P), (3.8) 

101 « (3.9) 

where is the infiaton field and V{(j)) is the scalar field potential. On defining the slow-roll 
parameters e and r]^ as |14U] 

H 

T]^ = e--^, (3.11) 
H(p 

the slow-roll conditions in Eqs. (13. 8p and (13. 9p translate into strong constraints for the 
slow-roll parameters: e, |?70| <ti 1, which actually become e, |?7^| < 10"^ in view of Eq. (13. 6 p 
for single-field inflation: 

nc- 1 = 27/^-6e, (3.12) 
and the observational requirements presented in Subsection | 



Multifield slow-roll models may also be characterized by a set of slow-roll parameters which 
generalize those in Eqs. (I3.10p and (13. lip |141j : 

e, = , (3.13) 



2 \V ^ 

r,, ^ ml^. (3.14) 

By writing the slow-roll parameters in terms of derivatives of the scalar potential, as in the 
latter two expressions, we realize that the slow-roll conditions require very flat potentials 
to be met. 



The level of non-gaussianity f^L in slow-roll inflationary models with canonical kinetic 
terms has been studied both for single-field |148] and for multiple-field inflation [221 1211[ 
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1226] . assuming ( series convergence and considering only the tree-level terms both in P,^ 
and in S^. What these works find is a strong dependence on the slow-roll parameters 
and rji] for instance, Ref. |148j gives us for single- field models: 

lfNL = e{l + f) + 2e-7]^, (3.15) 

where / is a function of the shape of the wavevectors triangle within the range < / < 5/6. 
Refs. [211 1138] show that in such a case the inclusion of loop corrections is unnecessary 
because the latter are so small compared to the tree-level terms. Thus, f^L in single-field 
models with canonical kinetic terms is slow-roll suppressed and, therefore, unobservably 
small. As regards the multifield models, Jnl was shown, first in the case of two- field 
infiation with separable potential |211] and later for multiple-field infiation with separable 
[22] and non-separable [226] potentials, to be a rather complex function of the slow-roll 
parameters and the scalar potential that in most of the cases ends up being slow-roll 
suppressed. Only for models with a sharply curved trajectory in field space might the Jnl 
be at most of order one, the only possible examples to date being the models of Refs. 
[31 [281 [Ml [37]. Again, such predictions are based on the assumptions that the ( series 
is convergent and that the tree-level terms are the leading ones, so they might be badly 
violated if loop corrections are considered. 

Following a parallel treatment to that in Ref. [211] . the level of non-gaussianity tnl is 
calculated in Ref. [188] for multifield slow-roll infiationary models with canonical kinetic 
terms, separable potential, and assuming convergence of the ( series and tree-level dom- 
inance. From reaching similar conclusions to those found for the /tvl case, the ttvl is 
slow-roll suppressed in most of the cases although it might be of order one if the trajectory 
in field space is sharply curved. Nevertheless, as we will shown in the next chapter, there 
may be a big enhancement in tnl if loop corrections are taken into account. 

Finally, it is worth mentioning that there are other classes of models where the levels of non- 
gaussianity /tvl and t^l are big enough to be observable. Some of these models correspond 
to general langrangians with non-canonical kinetic terms (/c-infiation [12], DBI infiation 
[195] , ghost infiation [11], etc.), where the sizeable levels of non-gaussianity have mostly a 
quantum origin, i.e. their origin relies on the quantum correlators of the field perturbations a 
few Hubble times after horizon exit. Non-gaussianity in has been studied in these models 
for the single- field case [HI 1187] and also for the multifield case [HI [HHl 1124] 1125] . A recent 
paper discusses the non-gaussianity in for these general models for single-field infiation 
[T5] . In contrast, there are some other models where the large non-gaussianities have their 
origin in the field dynamics at the end of infiation [261 1135] : nice examples of this proposal 
are studied for instance in Refs. [151] 1152] 1179] 1180] . However, since the infiationary models 
of the slow- roll variety with canonical kinetic terms are the simplest, the most popular, and 
the best studied so far although, in principle, the non-gaussianity statistical descriptors are 
too small to ever be observable, it is very interesting to consider the possibility of having 
an example of such models which does generate sizeable and observable values for f^i and 
Tnl- This appealing possibility will be the subject of the following sections and also the 
next chapter. 
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I SECTION sT 

A subclass of small-field slow-roll inflationary models 



According to the classification of inflationary models proposed in Ref . [57] , the small-field 
models are those of the form that would be expected as a result of spontaneous symmetry 
breaking, with a field initially near an unstable equilibrium point (usually taken to be at 
the origin) and rolling toward a stable minimum (0) ^ 0. Thus, inflation occurs when the 
field is small relative to its expectation value (j) ^ (0). Some interesting examples are the 
original models of new inflation [HI 1128] , modular inflation from string theory [55], natural 
inflation [51], and hilltop inflation [52] • As a result, the inflationary potential for small-field 
models may be taken as 

i 

where the subscript i here denotes the relevant quantities of the ith field, p is the same for 
all fields, and Aj and fii are the parameters describing the height and tilt of the potential 
of the ith field. 




While Ref. [2] studies the spectrum of ( for general values of the parameter p and an 
arbitrary number of fields, assuming ( series convergence and tree-level dominance, we will 
specialize to the p = 2 case for two fields and a: 



where we have traded the expressions 

A1+A2 for Vo, (3.18) 

^ for -Vo^, (3.19) 
/if 2mp 

and 

^ for -V^ottV- (3-20) 
/i2 2mp 

By doing this, and assuming that the first term in Eq. fl3.17p dominates, rj^ < and ?7cr < 
become the usual rj slow-roll parameters associated with the fields and a. 



We have chosen for simplicity the a = trajectory (see Fig. 13. ip since in that case the 
potential in Eq. f l3.17p reproduces for some number of e-folds the hybrid inflation scenario 
|129] where is the inflaton and a is the waterfall field. Non-gaussianity in such a model 
has been studied in Refs. [3l El [HSl [HSl |2071 [228]; in particular, Ref. [US] used a one- 
loop correction to conjecture that Jnl in this model would be sizeable only if ( was not 
generated during inflation, which turns out not to be a necessary requirement as we will 
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Figure 3.1: Our small-field slow-roll potential of Eq. f l3.17p with rj^,ricr < 0. The inflaton 
starts near the maximum and moves away from the origin following the a = trajectory 
depicted with the solid black line. (This figure has been taken from Ref. [5]). 

show later |l3l 1175] . Ref. [3], in contrast, works only at tree-level with the same potential 
as Eq. f l3.17p but relaxing the a = condition, finding that values for /at^ ~ are 
possible for a small set of initial conditions and assuming a saddle-point-like form for the 
potential (77^ < and ?7o- > 0). 



K^SECTION 3.5 



Constraints for having a reliable parameter space 



We will explore now the constraints that the model must satisfy before we calculate the 
level of non-gaussianity Jnl- Our guiding idea will be the consideration of the role that 
the tree-level terms and one-loop corrections in and have in the determination of the 
available parameter space. Only after calculating /tvl in Section |3T6] will we come back to 
the discussion of the consistency of the approach followed in the present section by studying 
the ( series convergence and the validity of the truncation at one loop level. 
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Tree-level or one-loop dominance: /] 



NL 



Since we are considering a slow-roll regime, the evolution of the background and a fields 
in such a case is given by the Klein-Gordon equation 

^ + 3H<i) + V^ = 0, (3.21) 

supplemented with the slow-roll condition in Eq. fl3.9p . This leads to 

cj>{N) = 0,exp(-iVr/^), (3.22) 

a{N) = a,exp{-N7]^), (3.23) 

so the potential above leads to the following derivatives of with respect to (j)^ and a^, for 
the a = trajectorylfl 



N, 



, N„„, 



-^exp[2N{ri^-ri^)] 
2vl 



,3^3 



exp[2N{r]^-7],)], 



(3.24) 
(3.25) 

0,(3.26) 



and so on. 

By means of the 5N formalism, we can make use of the above formulae to calculate the 
spectrum and the bispectrum of the curvature perturbation including the tree-level and 
the one-loop contributions when \ri„\ > {rj^l (see Appendix . This is the interesting case 
since, as will be shown in Section [3l6| it generates sizeable values for /atl. Following the 
results in Appendix |Al we will write down just the leading terms to the tree-level and one- 
loop contributions given in Eqs. ( lATT]) . ( lAl\ . (lA.lip and (lA.28p . 

2 



•jytree 



Tjtree 



2/A2 



El 

27r 



\n{kL) 



3M 



El 

2tx 



B. 



I— loop 

c 



exp[QN{r]^ - r]^)] [ ^ ] \n{kL)A7T 



(3.27) 



(3.28) 



(3.29) 



(3.30) 



Because of the exponential factors in Eqs. fl3.28p and (13.301) . it might be possible that the 
one- loop corrections dominate over and/or 5^. There are three posibilities in complete 
connection with the position of the field when the cosmologically relevant scales are 
exiting the horizon: 



^When calculating the A^-derivatives, we have considered that the final time corresponds to a slice of 
uniform energy density. This means, in the slow-roll approximation, that V is homogeneous. 
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Both Bi^ and P,^ are dominated by the one-loop corrections 

Comparing Eqs. f l3.27p with f l3.28p and Eqs. fl3.29p with fl3.30p we require in this case that 

4exp[4iV(|r/.|-|r/^|)] » — ^ , (3.31) 

^exp[6iV(|r/.|-|7^^|)] » ^ \ (3.32) 

in which case only the first inequality is required. Employing the definition for the tensor 
to scalar ratio r |137] : 

(3-33) 

1 /2 

Vj! being the amplitude of the spectrum for primordial gravitational waves, we can write 
such an inequality as 

« "-^i^Mmiv.l - \VM- (3.34) 
mpj 8 r]^ 

From now on we will name the parameter window described by Eq. fl3.34p as the low 0^ 
region since the latter represents a region of allowed values for 0^ limited by an upper 
bound. 



dominated by the one-loop correction and dominated by the tree-level 
term 



Comparing Eqs. ( I3.27P with ( I3.28P and Eqs. (I3.29P with fl3.30p we require in this case that 

^exp[4Ar(|r/.|-|r/^|)] « — ^ , (3.35) 

^exp[6iV(|r/<,|-|r/^|)] » — ^ , (3.36) 
''I' ^ [2^) 

which combines to give, employing the definition for the tensor to scalar ratio r introduced 
in Eq. (ICTll . 

'-^4^MmK\ - l^^l)] « f— y « '-^"^explQNiKl - . (3.37) 

From now on we will name the parameter window described by Eq. (13.371) as the interme- 
diate 0^ region since the latter represents a region of allowed values for 0^ limited by both 
an upper and a lower bound. 
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Both Bi^ and P,^ are dominated by the tree-level terms 

Comparing Eqs. f l3.27p with f l3.28p and Eqs. fl3.29p with fl3.30p we require in this case that 

^exp[m{\v^\-\v^\)] « ytW' ^^-^^^ 

^exp[6iV(|r/.|-|r/^|)] « ^ ^ (3.39) 

</.2 (27^) 

in which case only the second inequality is required. Employing the definition for the tensor 
to scalar ratio r introduced in Eq. fl3.33p . we can write such an inequality as 

1-)' » '-^$'Mm\v.\ - \VM . (3.40) 
mpj « 7]^ 

From now on we will name the parameter window described by Eq. (13.401) as the high 0^ 
region, since the latter represents a region of allowed values for 0^ limited by a lower bound. 



Spectrum normalisation condition 

The model must satisfy the COBE normalisation on the spectrum amplitude [35] 
considering that C is assumed to be generated during infiationj. There exist two possibilities 
discussed right below. 



( generated during inflation and dominated by the one-loop correction 



According to Eq. f l3.28p . and the tensor to scalar ratio r definition in Eq. fl3.33p . we have 
in this case 



which reduces to 



^exp[4iV(|r/.| - \^,\)] ln(A;L), (3.41) 



f— V= QVc4exp[4iV(|.^.|- 1.^,1)] ln(A;L), (3.42) 
\mpj \8J 7]l 



where Vt^ must be replaced by the observed value presented in Subsection 12.41 



^The scenario where ^ is assumed not to be generated during inflation will be presented in the next 
chapter. 
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( generated during inflation and dominated by the tree-level term 

According to Eq. fl3.27p . and the tensor to scalar ratio r definition in Eq. fl3.33p . we have 
in this case 



tree 



1 /i/*"' 



271 



V. 



which reduces to 

1 r 



1 f mp\^ rVt 

' ] (3-43) 



2 o • (3-44) 
mpj r]^8 

Notice that in such a situation, the value of the (p field when the cosmologically relevant 
scales are exiting the horizon depends exclusively on the tensor to scalar ratio r, once rj^ 
has been fixed by the spectral tilt constraint as we will see later. 



p. 5. 3 Spectral tilt constraint 

The combined 5-year WMAP + Type I Supernovae + Baryon Acoustic Oscillations data 
|119] constrain the value for the spectral tilt as 

n^- 1 = -0.040 ±0.014. (3.45) 

Here again we have two possibilities: is dominated either by the one-loop correction or 
by the tree-level term: 



dominated by the one-loop correction 

In this case the usual spectral index formula at tree-level |181] gets modified to account for 
the leading one-loop correction: 

n,-l = -Ae- 2ml ^;' ' + [ln(A:L)]"^ . (3.46) 

By making use of the derivatives in Eqs. fl3.24p . (13.250 . and (I3.26p . we have 

n^-l = -Ae + Ar]^+[\n{kL)]-\ (3.47) 

which implies that the observed value for is never reproduced in view of \n.{kL) ~ 0{1). 
Moreover, when calculating the running spectral index dnc^/dhik from Eq. (I3.47p . we 
obtain 

'^''^ - - [ln(A;L)]-' , (3.48) 



dhik 
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which rules out the possibihty that P,^ is dominated by the one-loop correction since the 
calculated dn^^/dlnk is far from the observationally allowed 2a range of values: —0.0728 < 
dn^/dlnk < 0.0087 [TTofl . 



dominated by the tree-level term 



Now the usual spectral index formula |181] applies: 



— 1 



(3.49) 



giving the following result once the derivatives in Eqs. f l3.24p . f l3.25p . and fl3.26p have been 
used: 

n^-l = -2e + 27]^ . (3.50) 

The efect of the e parameter may be discarded in the previous expression since, as often 
happens in small-field models [6], [32], e is negligible being much less than |?7o-|: 



■\t1<i>\ 













A < K\ 


\mpj 





(3.51) 



according to the prescription that the potential in Eq. f l3.17p is dominated by the constant 
term. Thus, by using the central value for — 1, we get 



1]^ = -0.020 



(3.52) 



fe.5.4- Amount of inflation 

Because of the characteristics of the inflationary potential in Eq. fl3.17p . inflation is eternal 
in this model. However, Ref. |157] introduced the multi-brid inflation idea of Refs. [1791 
1180] so that the potential in Eq. (I3.17P is achieved during inflation while a third fleld p 
acting as a waterfall fleld is stabilized in p = 0. During inflation p is heavy and it is trapped 
with a vacuum expectation value equal to zero, so neglecting it during inflation is a good 
approximation. The end of inflation comes when the effective mass of p becomes negative, 
which is possible to obtain if Vq in the potential of Eq. fl3.17p is replaced by 

K, = ^G(0,a)/ + ^(/-yy , (3.53) 

where 

G'(0,a) = (7l0=^ + ^72a^ (3.54) 



^We thank Eiichiro Komatsu for pointing out to us the dependence of n^^ and dn^/dlnfc on In(fcL). 
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S has some definite value and gi, g2, and A some coupling constants. The end of inflation 
actually happens on a hypersurface defined in general by |157[ 1180] 



In general the hypersurface in Eq. fl3.54p . defined by the end of inflation condition, is not 
a surface of uniform energy density. Because the 6N formalism requires the final slice to 
be of uniform energy densitjl^, we need to add a small correction term to the amount of 
expansion up to the surface where p is destabilised. In addition, the end of inflation is 
inhomogeneous, which generically leads to different predictions from those obtained during 
inflation for the spectral functions |^ I135[[T77] . In particular, large levels of non-gaussianity 
may be obtained by tunning the free parameters of the model. Specifically, by making 
gi/92 ^ I5 large values for /tvl and t^l are obtained due to the end of inflation mechanism 
rather than due to the dynamics during slow- roll inflation [U [981 1157] . 

Ref. [38] chose instead the case gf/g2 = V4>/VfT such that the surface where p is destabilised 
corresponds to a surface of uniform energy densitj0. In this case all of the spectral functions 
are the same as those calculated in this thesis (see Refs. HSl 1175] ). which in turn are 
valid at the final hypersurface of uniform energy density during slow-roll inflation. Thus, 
we have a definite mechanism to end inflation which, nevertheless, leaves intact the non- 
gaussianity generated during inflation. 

It is well known that the number of e-folds of expansion from the time the cosmological 
scales exit the horizon to the end of inflation is presumably around but less than 62 
[MH [220] • The slow-roll evolution of the (f) field in Eq. IKT2^ tells us that such an 
amount of inflation is given by 



?70 \ 0* y ~ 

where (pend is the value of the field at the end of inflation. Such a value depends noticeably 
on the coupling constants in Eq. f l3.53p . We will in this chapter not concentrate on the 
allowed parameter window for gi, g2, and A. Instead, we will give an upper bound on 
during inflation, for the 77^ < case, consistent with the potential in Eq. f l3.17p and the 
end of inflation mechanism described above. 

Keeping in mind the results of Ref. [13] which say that the ultraviolet cutoff in cosmological 
perturbation theory could be a few orders of magnitude bigger than mp, we will tune 
the coupling constants in Eq. (13. 53^ so that inflation for rj^ < comes to an end when 
|?7(^|</)^/2mp ~ 10"^. This allows us to be on the safe side (avoiding large modifications 
to the potential coming from ultraviolet cutoff-suppressed non-renormalisable terms, and 
keeping the potential dominated by the constant Vq term). Coming back to Eq. (I3.56p . we 



^See for instance Ref. [35|. 

'^Ref. [55] studies as well the case where gf = g|, but we are not going to consider it here. 



G(0,a) = S2. 



(3.55) 




(3.56) 
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get then 



which leads to 



N 



In 



2 X 10^2 y/2 



< 62, 



> 



nip 



2 X 10- 



2\ 1/2 



exp(-62|r/<^|) . 



(3.57) 



(3.58) 



.SECTION 3.6. 



Non-Gaussianity: /^vl 



In this section we will calculate the level of non-gaussianity represented in the parameter 
Inl by taking into account the constraints presented in Subsections I3.5.2[ I3.5.3[ and 13. 5. 4t 
and the different (j)^ regions discussed in Subsection 13.5.11 



The low 0^ region 

This case is of no observational interest because P(; dominated by the one-loop correction 
is already ruled out by the observed spectral index and its running as shown in Subsection 
13.5.31 In addition, the generated non-gaussianity is so big that it causes violation of the 
observational constraint f^L > —9: 

„ j^l~loop 

5^- = 4,,E.4,^.-o^^. = -In>n(«r"= ~ -2 X 10'. (3.59) 

according to the expressions in Eqs. fl2.35p . f l2.3ip . f l3.28p . and fl3.30p . 

We want to remark that, although it is of no observational interest, this case represents the 
first example of large non- gaussianity in the bispectrum of C, for a slow-roll model of 
inflation with canonical kinetic terms. It is funny to realize that the model in this case is 
additionally ruled out because the observational constraint on /jvl is violated by an excess 
and not by a shortfall as is currently thought [221 EHl [Ml EIB 1226] . 
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The intermediate 0^ region 

The level of non-gaussianity, according to the expressions in Eqs. fl2.35p . f l2.3ip . fl3.27p . and 
fl3.3Up . is in this case given by 

r^l-loop 3 / tr \ 2 



^fML = = ^exp[6iV(|,7.|-k<^l)]ff^) HkL) 

^ eMQN{\va\ - M] (^) ' ^ 



= r?^exp[6Ar(|??,|-|r?^|)]P<^ln(fcL), (3.60) 

^ tInl ~ -2.457 X 10~^|?7„|^exp[300 ln(5.657 X (|r?^| -0.020)] , (3.61) 

o 

where in the last line we have used expressions from Eqs. fl3.44p . fl3.52p . and f l3.57p . 

Now, by implementing the spectral tilt constraint in Eq. fl3.52p in the spectrum normal- 
isation constraint in Eq. f l3.44p and the amount of inflation constraint in Eq. f l3.58p . we 
conclude that the tensor to scalar ratio r is bounded from below: r > 2.680 x 10~^. 



In the plot r vs |?7o-| in figure [3l2| we show lines of constant f^L corresponding to the values 
/jvL = ~5, —10, —15. We also show the high and intermediate 0^, regions in agreement with 
the constraint in Eq. (I3.37P : 

8 r]^ \mpj 8 7]^ 

8.139x10^ < |r/^pexp[3001n(5.657xl0-V-i/2) (^|^^| _ o.020)] < 8.210x10^^ (3.62) 

As is evident from the plot, the WMAP (and also PLANCK) observationally allowed 2a 
range of values for negative /nl, —9 < /nl, is completely inside the intermediate 0^, region 
as required. More negative values for f^i, up to /nl = —20.647 are consistent within 
our framework for the intermediate 0* region, but they are ruled out from observation. 
Nevertheless, like for the low (j)^ region studied above, it is interesting to see a slow- roll 
inflationary model with canonical kinetic terms where the observational restriction on /jvl 
may be violated by an excess and not by a shortfall. So we conclude that if is dominated 
by the one-loop correction butPc^ is dominated by the tree-level term, sizeable non-gaussianity 
is generated even if C, is generated during inflation. We also conclude, from looking at the 
small values that the tensor to scalar ratio r takes in figure 13.21 compared with the present 
technological bound r > 10~^ [65], that for non-gaussianity to be observable in this model, 
primordial gravitational waves must be undetectable. 

Notice that in oder to get positive values for fNL-, which is observationally more interesting 
in view of the results presented in Subsection 12. 4[ ri„ should be positive according to Eq. 
fl3.60p . However, being rj^ negative in order to reproduce the observed spectral tilt, the 
argument of the exponential in Eq. (I3.60p would be negative, making the fNL obtained too 
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0.092 0.094 0.096 0.098 0.100 



\ria\ 

Figure 3.2: Contours of Jnl in the r vs |?7o-| plot. The intermediate (high) 0^ region 
corresponds to the shaded (white) region. The WMAP (and also PLANCK) observationally 
allowed 2a range of values for negative /atl, — 9 < /atl, is completely inside the intermediate 
0* region. Notice that the boundary line between the high and the intermediate 0* regions 
matches almost exactly the fNL = —1.667 line. 

small to be observationally interesting. As regards the general case, in view of the previous 
reason being model dependent, we may only say that in order to get Jnl positive when 5^ 
is dominated by the one-loop corrections, Bi^ should be positive (based on the definition 
of fNL in Eq. f l2.3ip ) which means that the maximum between N^a and iV^^ should be 
positive in view of Eq. f lA.12|) . 

Finally we want to point out that, by reducing our model to the single-field case, the 
consistency relation between /tvl and presented in Ref. |16]: /nl ~ ^{^c ~ 1) is not 
violated since in that case is never dominated by the one-loop corrections for slow-roll 
inflation as demonstrated in Ref. [1380. Thus, the level of non-gaussianity /nl for our 
model reduced to the single-field case is described by the high 0^ region as shown below. 



^We thank Eiichiro Koniatsu for questioning us about this issue. 
^We thank Fihppo Vernizzi for questioning us about this issue. 
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The high 0^ region 

This case is of no observational interest because, according to the expressions in Eqs. 
fl2.35p . fl2.3ip . f l3.27p . fl3.29p . and fl3.52p . the non-gaussianity generated is too small to be 
observable: 

tInl = = -V<t> = 0.020 , (3.63) 

47r4^(Pf-)2 

in agreement with the consistency relation of Ref. [16] for our model reduced to the single- 
field case, and with the general expectations of Refs. [221 1148[ I19H I21H 1226] for slow-roll 
inflationary models with canonical kinetic terms where only the tree-level contributions are 
considered. 



SECTION 3.7 



Convergence of the ( series and perturbative regime 



In Sections 13.51 and 13.61 we have worked up to the one-loop diagrams in order to constrain 
the parameter space and find the level of non-gaussianity /tvl- It is time then to address 
the issue of the ( series convergence and justify the existence of a perturbative regime so 
that the truncation of the series up to the one-loop order, for the model we have considered, 
is valid. A way to do that is by rederiving the ( series in terms of and 5(T^ by equating 
the unperturbed scalar potential to the perturbed one at the final time t; this of course is 
valid in view of the first slow-roll condition in Eq. ( 13. Sp and the final slice being one of 
uniform energy density: 

\^ Zi if L p Zi I lb p 

= vAl + ^r/<^^^^^±|^exp[-2(iV + 5N)r^^] + ir/.^^^±|^ exp[-2(iV + 5N)r^A . 

(3.64) 

From the previous expression it follows that 

?7<^0^exp[-2iV?7^] + r]^alexp[-2Nr]^] 

= 7]^{(f), + 50,)2 exp[-2(iV + 6N)r]^] + 7]^{a, + Sa^f exp[-2(iV + SN)r],] , (3.65) 
which is easier to handle in terms of variables x and y defined as 

X = ^, (3.66) 



y = 



fi + — 



exp[2iV(?70 - r]^)] 



1/2 



(3.67) 



Convergence of the ^ series and perturbative regime 



41 



Thus, the exponentials factors contaning N (but not 6N) are completely absorbed in y, 
and the expression in Eq. fl3.65p looks as follows: 



1 + 



= (1 + xf exp[-25iVr/^] + "^y^ eM-25Nr^„ 



(3.68) 



If we were able to solve for SN in Eq. ( I3.68P in terms of r]^, r^o-, x, and y (after making 
cr^ = 0), we could Taylor-expand around x = and y = reproducing the x-dependent 
part of Eq. fl3.3p . This would be really good because the Taylor expansion would look so 
clean, in the sense that all the concerning exponential factors contaning which appear 
explicitely in Eq. (13. 3p would already be absorbed in y, that the issue of truncating at some 
specific order in 50* and ^a* would be simply justified by requiring |x| ^ 1 and |?/| ^ 1. 
Nevertheless, as is seen in Eq. (13.681] , it is impossible to solve for 6N in terms of 77^, rj^, x, 
and y unless we make a Taylor expansion of the exponential functions aroud 6N = 0: 



[(1 + 



X 



1] + 



-6N 



-2v,{l + x f - 2%' 



1 + ^ 



6N' [2r/5(l 



X) 



(3.69) 



Notice that the Taylor expansion of the exponential functions is always convergent whatever 
the arguments of the exponentials are. Moreover, if the Taylor expansion derived from a 
function f{x) converges, it converges precisely to f{x) [200] . Thus, the expression in Eq. 
(I3.69P is actually the same as the expression in Eq. (I3.68p . 



Now, solving for 6N in terms of 77^, r]cr, x, and y, although possible in view of the expression 
in Eq. (I3.69p . is not an easy business. That is why we will truncate the series in Eq. (I3.69p 
up to second order in 5N and solve the resultant quadratic equatioij^. Notice that, since 
( = 6N — (SN) and ( ~ 10~^, we may truncate the series in Eq. (I3.69P up to whatever 
order we wish and still reproduce ( to high accuracy. Thus, the solution for the quadratic 
equation coming from the series in Eq. (I3.69P after truncation at second order is: 



SN 




1] + 



[r,^ [{l+xf + y']Y\ 



1 



1 + ^ 



[(1 



X] 



(3.70) 



"'^'^The truncation up to second order in SN has been chosen in order to have complete consistency with 
the order of the variables x and y in Eq. p.68p . 
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If in addition we make Taylor expansions of the square root and the factor in the third hne 
of the previous expression around x = and y = 0: 



X] 



x] 



1 



+ 

2 2ri^ 



1 _ m 



1 + ^ 



1] + 



1 + ^ 



1/2 



[r]^ [(1 + xY + y^]} ' = — [1 - 2x + 3x 



V<t> 



...] 



introducing them into Eq. (]3.70p . we end up with the following power series for 6N: 



6N 



X 



X 



2^ 



(3.73) 



where the ± symbol is changed to the — sign so that 6N remains a perturbation, and the 
trajectory cr = is chosen. Coming back to the variables 6(f)^ and Sa^, we see that Eq. 
(13.731) reproduces the x-dependent part of Eq. (13. 3p in view of Eqs. (I3.24p and (I3.25P up 
to second order in Sch^, and 5cr^,: 



6N 



Eq. (I3.73p . although reliable only up to second order, tells us that the expected behaviour of 
6N in terms of r]^, rjcr, x, and y is indeed obtained. Moreover, from our previous discussion 
we know that 6N can be exactly written in terms of a series of x and y withouth the 
explicit appearance of the concerning exponential factors containing A^. This is indeed 
partially confirmed up to third order when introducing Eqs. (I3.24p . (I3.25p . and (I3.26P into 
the x-dependent part of Eq. (13. 3p : 



6N 



X 



X 

Y 



2 2/ 



X'^ 

y 



V4' 2 
—^xy 



(3.75) 



The bottom line of this discussion is that we have been able to identify two quantities that 
determine the truncation of the series up to some specific order. These two quantities are 
X and y which we could identify as the "coupling constants" of the theory in the context of 
Quantum Field Theory. By making |a;| <C 1 and \y\ <C 1 we can see from Eq. (I3.75P that all 
the terms higher than second order in x and y are subleading compared to the second-order 
ones. As regards the first-order terms compared to the second-order ones, we see that the 
latter are not necessarily subleading compared to the former because of the non-existence 
of the first-order y term and in view of \y/x\ < 1600 from Eqs. (I3.66P and (I3.67P and the 
values for 77,^, -q^ and considered in Sections 13. 5[ 13.61 and 14.61 In the language of the 
Feynman-like diagrams [39], truncating 6N in Eq. (I3.75P up to second order in x and y 



(3.71) 
(3.72) 
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means considering only the leading diagrams at tree level and one loop which is what we 
have done in Sections 13. 5[ 13.61 In fact, |x| <C 1 and \y\ <C 1 mean that 



\x\ 




(3.77) 

which are well satisfied for the cases when is dominated by the tree-level term (see 
Subsection EXl- Eq. (K^ and Subsection EXH- Eq. IK^ ): 



2ti 




\Vl''^ ^ 10-^ , (3.78) 



||(§)'^exp[2iV(r/^-r/.)]| « |^ exp[-2iV(|r/.| - |r/,|)] j'^' <2. (3.79) 



By explicitly calculating the two-loop and three- loop diagrams for and P^, and employing 
the results of Ref. |lUUj . we have checked that the conditions |x| ^ 1 and |?/| ^ 1 efectively 
make these diagrams subleading compared to the leading ones at one-loop level. 

Finally, we will discuss the convergence of the C series in view of Eqs. (13. 70 p . (13. 71 p . and 
(13.721) . We first note that the series in Eq. (13.711) is always convergent. As regards the 
series in Eq. (13. 72 p . it will not be convergent at all while the function 

1 1 



+ + ^ (l + x)2 + p2x2 ' 
with 

3/2 



(3.80) 



exp[N{r]^-r],)], (3.81) 



does not satisfy the following necessary condition |200] : for the Taylor series around x = 
of a function /(x) to be convergent, it is necessary that the extension f{z) to the complex 
plane of f{x) is continous in a neighbourhood of z = 0. If this is the case, and the Taylor 
series of f{z) is indeed convergent, the convergence circle must either match or be inside 
the aforementioned neighbourhood. Of course, this is not a sufficient condition, but at least 
gives us a constraint on the possible values that x may take. 

Applying this condition to the expression in Eq. (I3.80p . we see that the extension of this 
function to the complex plane has poles for (1 + z)"^ = —B^z^ which leads to 



_ ±'iB - 1 

'~ P2 + l ■ 

Therefore, the extension to the complex plane of Eq. (I3.80p is continous for 



(3.82) 



kl<|^ = l' (3-83) 
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so the necessary condition for the convergence of the series in Eq. (13 .72^ . and therefore 
for the convergence of the series in Eq. (13.701) which is what we are interested in, is 
given by \x\ < 1. Thus, such a necessary condition for the convergence of the ( series is 
automatically satisfied once we choose \x\ <^ 1, as we have seen above it is required for 
working in a perturbative regime. 

I SECTION 3.8^ 

Conclusions 



Is it reasonable to study the primordial curvature perturbation ( by identifying it with 
a truncated 6N series expansion? Is it actually possible to cut out with confidence such 
a series at some specific order? Is it true that all the slow-roll inflationary models with 
canonical kinetic terms produce primordial non-gaussianity supressed by the slow-roll pa- 
rameters? Are the loop corrections in cosmological perturbation theory always smaller than 
the tree-level terms? We have addressed these questions in this chapter, answering all of 
them by paying particular attention to a special slow-roll inflationary model. The ( series 
expansion is indeed a powerful tool to study the statistical descriptors of (; nevertheless, 
we should seek for the convergence radius in order not to obtain results that actually have 
nothing to do with (. We may cut out the series but, to be completely sure about the 
precision of our approximations, we have to study the conditions for the existence of a per- 
turbative regime. Non-gaussianity in slow-roll inflationary models with canonical kinetic 
terms is not always suppressed by the slow-roll parameters; we have seen this at tree-level 
for /jvi in Refs. [31 EZ] , and considering loop corrections for f^^ in the prersent chapter 
and in Refs. |l3l |38]. Particulary, we shown in this chapter that it is possible to attain 
very high, including observable, values for the level of non-gaussianity /nl associated with 
the bispectrum 5^ of the primordial curvature perturbation (, in a subclass of small-field 
slow-roll models of inflation with canonical kinetic terms. Such a result was obtained by 
taking care of loop corrections both in the spectrum and the bispectrum and as- 
suming that the latter can dominate over the former; of course, this possibility is model 
dependent. More precisely, we can say that if 5^ is dominated by the one- loop correction 
but is dominated by the tree-level term, sizeable non-gaussianity is generated even if ( 
is generated during inflation. What is interesting is that this kind of particular models are 
populary known to predict too small values for the level of non-gaussianity /nl, as small 
as the slow-roll parameters. Finally, as far as we have investigated, the loop corrections in 
cosmological perturbation theory are not always smaller than the tree-level terms; in fact, 
when they become the leading contributions, a surprising phenomenology appears in front 
of our eyes. 



Chapter 




ON THE ISSUE OF THE ( SERIES 
CONVERGENCE AND LOOP CORRECTIONS IN 
THE GENERATION OF OBSERVABLE 
PRIMORDIAL NON- GAUSSIANITY IN 
SLOW-ROLL INFLATION: THE TRISPECTRUM 



SECTION 4.1 

Introduction 



The primordial curvature perturbation ( [56 |I140[ [155(1220] . and its SN expansioE0 [52lll39[ 
I142[ I18H I182[ 1202] . was the subject of study in a previous chapter (see also [13]). We were 
interested in how well the convergence of the ( series was understood, and if the traditional 
arguments to cut out the ( series at second order [142] 1228] . keeping only the tree-level 
terms to study the statistical descriptors of C [221 EH SHI [ffl I22SI [223 [22Z], were 
reliablell. We argued that a previous study of the ( series convergence, the viability of a 
perturbative regime, and the relative weight of the loop contributions against the tree- level 



^By ^^6N expansion" we mean approximating 5N by a power series expansion in the initial conditions. By 
'^SN formula" we mean the statement that to lowest order in spatial gradients ( = 6N. These conventions 
have been and will be used throughout the text. 

^We follow the terminology of Ref. [31] to identify the tree-level terms and the loop contributions in a 
diagrammatic approach. The associated diagrams are called Feynman-like diagrams. 
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terms, were completely necessary and in some cases surprising. For instance, the levels of 
non-gaussianity /nl and tnl in the bispectrum i?^ and trispectrum of ( respectively, 
for slow- roll inflationary models with canonical kinetic terms |137[ 1 1401 1141] . are usually 
thought to be of order 0(ei, r/i) [221 12111 1226^ and 0(r) |l88l Il92j3respectively, were 



and rji are the slow-roll parameters with e^, |?7j| <C 1 |141] and r is the tensor to scalar ratio 
|137] with r < 0.22 at the 95% confidence level [119] . However, in order to reach such 
a conclusion, generic models were used where the loop contributions are comparatively 
suppressed and, therefore, the truncated 6N expansion may be used. Of course exceptions 
may occur, and in those cases it is crucial to check up to what order the truncated 6N 
expansion may be used, and which loop contributions are larger than the tree-level terms. In 
any of these cases, general models or exceptions, the question regarding the representation 
of ( by the SN expansion is a matter to discuss. 

Refs. [H [37] show that large, and observable, non-gaussianity in is indeed possible for 
certain classes of slow-roll models with canonical kinetic terms and special trajectories in 
field space, relying only on the tree-level terms. Ref. [3B] does the same for and but 
this time arguing that the loop corrections are always suppressed against the tree-level terms 
if the quantum fluctuations of the scalar fields do not overwhelm the classical evolution. 
Nonetheless, although the resultant phenomenology from papers in Refs. [5l [371 [38] is very 
interesting, the classicality argument used in Ref. [38] is very conservatively stated leading 
to too strong conclusions as we will argue later in this chapter. More research remains 
to be done to understand the role of the quantum diffusion and, being this beyond the 
scope of the present chapter, we will leave the discussion for a future research project. 
We addressed the ( series convergence and the existence of a perturbative regime in the 
previous chapter, showing how important the requirements to guarantee those conditions 
are. Moreover, we showed that for a subclass of small-field slow-roll inflationary models 
with canonical kinetic terms, the one-loop correction to 5^ might be much larger than 
the tree-level terms, giving as a result large, and observable, non-gaussianity parameterised 
by Jnl- The present chapter extends the analysis presented in the previous one to 
showing, for the first time, that large and observable non-gaussianity parameterised by ttvl 
is possible in slow-roll inflationary models with canonical kinetic terms due to loop effects, 
in total contrast with the usual belief based on the results of Refs. |1881 1192] . In order to 
properly identify the non-gaussianity levels found in previous chapter and in the present 
one with those that are constrained by observation, we comment on the probability that 
an observer in an ensemble of realizations of the density field in our scenario sees a non- 
gaussian distribution. As we will show such a probability is non-negligible for the concave 
downward potential, making indeed the observation of the non-gaussianity studied in this 
chapter quite possible. 

The layout of the chapter is the following: in Section [^^ we make some aditional comments 
about ot the slow-roll inflationary model that exhibits large levels of non-gaussianity when 



3 See however Refs. [SJ|37]. 
-^See however Refs. [raSlIgS]. 
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loop corrections are considered. This model was described in more detail in Section 13.41 
In Section 14.31 we study the impact of the quantum Auctions of the scalar fields on their 
classical evolution. As a result we argue how the loop suppression proof given in Ref. [38] 
does not apply to our model. Section 14.41 studies the probability of realizing the scenario 
proposed in this thesis for a typical observer. Section 14.51 is devoted to the reduction of 
the available parameter window by taking into account some restrictions of methodological 
and physical nature. The level of non-gaussianity tnl in the trispectrum is calculated 
in Section 14.61 for models where ( is, or is not, generated during inflation; a comparison 
with the current literature and the results found in the previous chapter for f^i is done. In 
Section B. 7. 2l the level of non-gaussianity f^L in the bispectrum is calculated for models 
where ( is not generated during inflation. Finally, Section 14.81 presents the conclusions. 

— ^« j;F,r;TTmv d.-?. 

A quadratic two-field slow-roll model of inflation 



We will give in this Section some relevants remarks about the inflationary potential studied 
in Section 13.41 and given by 



where and a are the inflaton fields and mp is the reduced Planck mass. By assuming that 
the first term in Eq. (14.11) dominates, r]^ and rj^ become the usual rj slow-roll parameters 
associated with the fields and a. 

We have chosen the a = trajectory since this case is the easiest to work from the point 
of view of the analytical calculations and because it gives the most interesting results. In 
addition, such a trajectory for the potential in Eq. (14. ip reproduces for some number of 
e-folds (for ?70,?7o- > 0) the hybrid inflation scenario |129] where is the inflaton and a is 
the waterfall field. We will analyze in Section I4l4l the probability for an observer to live in a 
region where a = for the concave downward potential. Non-gaussianity in the bispectrum 
Bt; of C for this kind of model has been studied in Refs. EH ESI SSI EQ M ESB] ; 
in particular, Ref. [13] shows that the one-loop correction dominates over the tree-level 
terms if rj^^rjfj < and |?7o-| > |?70|, generating in this way large values for Jnl even if ( 
is generated during inflation. Refs. [3 |37], in contrast, work only at tree-level with the 
same potential as Eq. (14. ip but relaxing the a = condition, finding that large values for 
/tvl are possible for a small set of initial conditions. Ref. [5B] improves the analysis in 
Ref. [3Z], this time taking into account also the trispectrum of ( and the role of the 
loop corrections. According to that reference, large values for tnl are also possible for a 
small set of initial conditions if the tree-level terms dominate over the loop corrections. 
Moreover, it is claimed that loop corrections for this model are always suppressed against 




(4.1) 
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the tree-level terms if the quantum fluctuations of the fields are subdominant against their 
classical evolution. The opposite case seems to happen for some narrow range of initial 
conditions including a* = 0, which is the case studied in this chapter. As we will argue 

is expressed in a very conservative 



in Section 14. 3[ the classicality condition in Ref. [38] 
way, leading to too strong and non-general conclusions. Dominance of loop corrections is 
then safe from the classical vs quantum condition, allowing the interesting large levels of 
non-gaussianity discussed in Chapter [3] and in the present one. 



Following the results in Appendix |X1 we will write down the leading terms to the spectrum, 
bispectrum, and trispectrum of the primordial curvature perturbation ( including the tree- 
level and one-loop contributions given in Eqs. f lXTTj) . (E2SD, (IS2SD, dSSHD, flX29|) and 
(1X301) : 
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(4.7) 



where L is the infrared cutoff chosen so that the quantities are calculated in a minimal 
box [271 1138] . Except when considering low CMB multipoles, the box size should be set at 
L ^ Hq [1141 1121] . giving \n{kL) ~ 0{1) for relevant cosmological scales. 



The important factor in the loop corrections is the exponential. This exponential function 
is directly related to the quadratic form of the potential with a leading constant term. It 
will give a large contribution if rj^ > rjcr- In Chapter [3l we chose the concave downward 
potential in order to satisfy the spectral tilt constraint, which makes i]^ < 0, while keeping 
\Vu\ > \v<l>\- III this chapter we will consider the same case. 
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- I SECTION IT 

Classicality 



Ref. [38] argues in Appendices A and B how, by imposing the requirement that the quantum 
fluctuations of the fields around their background values do not overwhelm the respective 
classical evolutions, the loop corrections to P^, -B^, and are suppressed against the tree- 
level terms. The proof relies on the fact that, if the classicality condition is satisfied, the 
second-order terms in the SN expansion in Eq. fl3.3p . are subleading against the first-order 
terms. This in turn implies 

pi— loop 

4^7 <1, (4.8) 

pl—loop 

« 1 ' (4-9) 

(4-10) 

as explicitly stated in Eqs. (A.16-A.19) of Ref. [38]. In addition, under the same as- 
sumptions, higher order corrections in the spectral functions P^, P^, and are always 
subleading against the one-loop corrections and, therefore, subleading against the tree-level 
terms. This conclusion is obtained if the 5N expansion may be truncated at fourth order. 
However, what is the classicality condition employed in Ref. [38] ? 

Assuming slow-roll evolution for each field 0j, which is valid only if the quantum fluctuation 
5(f)i is by far smaller than the classical evolution A0j, the classical change in the ipi field 
during a Hubble time around horizon exit is 

A^.W=»-Jgg. (4.11) 

where Vi denotes the derivative of the potential with respect to the i-th field. Comparing 
the latter expression with the quantum fluctuation 

(4.12) 

and requiring that A0j is much larger than we get 



For our quadratic two-field slow-roll model of inflation, where the slow-roll evolution is 
given by Eqs. f l3.22p -( l3l23|) . the previous expression translates into 



P* 



(4.14) 
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which is the one given in Eq. (A.l) of Ref. [38]. Such a condition is equivalent to 



6(f)i 



V6 



(4.15) 



which is the one given in Eq. (A. 2) of Ref. [38]. Thus, under this condition, the trajectory 
a = studied in this thesis seems not to be well described by the slow-roll approximation 
and, therefore, the obtained results based in the 6N formalism would not be reliable. 



This classicality argument given in Ref. [3B] is too conservatively stated. To see why it is like 
that, we may reason in the following way for general inflationary models: for any point along 
the background classical trajectory in field space it is possible to rotate the field axes so that, 
instantaneously, there is an infiaton (or 'adiabatic') field that points along the trajectory 
and some light 'entropy' fields which point in orthogonal directions [77]. The quantum 
fluctuations for the entropy fields are non-vanishing but the classical evolution for each of 
these fields is zero. Since the condition in Eq. fl4.13p is not formulated in any particular field 
parameterisation, we may argue that for any multi-field inflationary model the application 
of this condition would lead to a background inflationary trajectory dominated by the 
quantum evolution. Thus, slow-roll conditions would always be impossible to apply. A 
more general classicality condition should still be the one in Eq. fl4.13p but only applied 
to the adiabatic fleld and not to the entropy fields. In that respect the a = trajectory 
studied in this thesis is safe from large quantum fiuctuations since the classicality condition 
in Eq. fl4.13p applied only to the field is extremely well satisfied as long as <^ 1, which 
is actually the case for single field slow-roll infiation |174] . Indeed, a much better way of 
stating the classicality condition is the foUowingjfl if the infiationary trajectory must be 
dominated by the classical motion of the fields, then the perturbation in the amount of 
infiation, due to the quantum fiuctuations of the fields, must be negligible: 

(5iV < 1 . (4.16) 

By virtue of the 6N formalism, this expression is simply satisfied if the free parameters 
of the infiationary model under consideration are chosen so that the COBE normalisation 
(V^^"^ ^ 5 X 10~^ [35]) is satisfied, which is always the case. Nevertheless, we understand 
that the role of quantum diffusion is of great importance (see for instance Refs. [SH 1168] ). 
and a dedicated study of this issue is left for a future research project. 



^We acknowledge Misao Sasaki for pointing out to us this idea. 
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— I section"1^ 
Probability 



The main purpose of this chapter is to identify regions in the parameter space with high 
levels of primordial non-gaussianity. Then, we proceed to compare the obtained non- 
gaussianity with observation. In order to do the latter, we first need to realize what the 
probability is for a typical observer to live in a universe where the inflationary trajectory is 
the one studied in this thesis: a = 0. This is particularly relevant for the concave downward 
potential where the background trajectory a = is unstable. 

In the context of quantum cosmology, the probability of quantum creation of a closed 



which means that the universe can be created if V is not too much smaller than the Planck 
density. Thus, for our concave downward potential, having chosen the field contributions 
to V in Eq. (I3.17P to be negligible is good because it increases the probability. In addition, 
within a set of initial conditions for 0, the most probable initial condition for a is a = 0. 
The (j) = trajectory is also highly probable but, since we are assuming |?7o-| > Ir/^l, the 
a = trajectory is more probable. This of course implies that the levels of non-gaussianity 
obtained in this thesis may be observable. 

^SECTION 4.5' 

Reducing the available parameter window 



The analysis of the observed spectral index and the ( series convergence was given in 
Subsection 13.5.21 and in Section 13. 7[ respectively. Regarding the existence of a perturbative 
regime, we have to add to the discussion in Section [3171 that, by cutting out the series in 
Eq. fl3.75p at second-order, just one Feynman-like diagram per spectral function of ( is 
necessary to study the loop corrections to these spectral function^. That is why in the 



^This is assuming that the Feynman-hke diagrams containing n-point correlators of the field pertur- 
bations with n > 3 are subdominat against the diagrams containing only two-point correlators. For the 
trispectrum this does not happen when the tree- level terms dominate over the loop corrections |188[ 1192] . 
However, for the cases considered in this chapter, when the loop corrections in the trispectrum dominate 



universe is proportional to 12121 HH] 




(4.17) 
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previous chapter there was just one leading diagram for the one-loop correction to (Fig. 
IA.2b ). as well as one leading diagram for the one-loop correction to (Fig. IA.4b ). When 
applied to T^, this analysis shows that the only diagrams to consider are the one in Fig. 
IA.5b for the tree-level terms, and the one in Fig. IA.5b for the loop corrections. Such 
diagrams lead to the expressions in Eqs. f0r29D and fOOOD for T*''^^ and T^-^°°p. 

In the following, we will give the relevant information for when ( is generated during 
inflation, and for both 5^ and when ( is not generated during inflation. 



^^^1 Tree-level or one-loop dominance: tnl 

The exponencial factors in Eqs. (14. 3 p and (14. 7p open up the possibility that the loop 
corrections dominate over "P^ and/or T^. There are three posibilities: 



Both and P^ are dominated by the one-loop corrections 

Comparing Eqs. ( 14. 2p with ( 14. 3p and Eqs. (14. 6 p with ( 14. 7p we requiere in this case that 

4^exp[8iV(r/^-r/.)] » ^ ^ , (4.19) 

in which case only the first inequality is required. Employing the definition for the tensor 
to scalar ratio r introduced in Eq. (I3.33p . we can write such inequahty as 

f — y « ^4exp[4iV(r/, - V.)] ■ (4.20) 
\mpj 8 r]j 

From now on we will name the parameter window described by Eq. (I4.20p as the low 0^ 
T-regiorQ, since the latter represents a region of allowed values for 0^ limited by an upper 
bound. 



over the tree- level terms generating in turn large values for tn l , the diagrams containing n-point correlators 
of the field perturbations with n > 3 are expected to be subdominat because of their dependence on the 
slow-roll parameters |100j . 

''The T in T-region is introduced in order to differentiate explicitly between these regions and those 
found in the subsection 13 . 5 . II for B^. 
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dominated by the one-loop correction and dominated by the tree-level 
term 



Comparing Eqs. (14.21) with (14. 3 p and Eqs. (14. 6 p with (14. 7p we requiere in this case that 

^exp[4iV(r/^-r^.)] « ' ^^'^^^ 

Aexp[8N{v^-V^)] » — (4.22) 

which combine to give, employing the definition for the tensor to scalar ratio r introduced 
in Eq. (K33\\ . 

^1 exp[4iV(ry, - r^^)] « (^-^^ ' « exp[8iV(ry, - r^.)] . (4.23) 

From now on we will name the parameter window described by Eq. (14.231) as the interme- 
diate 0^ T-region, since the latter represents a region of allowed values for 0* limited by 
both an upper and a lower bound. 



Both and are dominated by the tree-level terms 

Comparing Eqs. (14.21) with (14. 3 P and Eqs. (14. 6 P with (14. 7p we requiere in this case that 

^exp[4iV(r/^-r/.)] « ' ^^'^^^ 

4!|exp[8iV(r/^-r/.)] « ^ , (4.25) 

'^'f' 02 ( 2^ j 

in which case only the second inequality is required. Employing the definition for the tensor 
to scalar ratio r introduced in Eq. (I3.33p . we can write such an inequality as 

» '-^^^Mmv, - V.)] . (4.26) 

From now on we will name the parameter window described by Eq. (14.260 as the high 0* 
T-region, since the latter represents a region of allowed values for limited by a lower 
bound. 



Ik. 5 .2 The normalisation of the spectrum 



Either ( is or is not generated during infiation, we must satisfy the appropriate spectrum 
normalisation condition. There exists four possibilities; however, it was shown in Subsection 
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I3.5.2l that the case where ( is generated during inflation and is dominated by the one-loop 
correction, is no of observational interest since it is impossible to reproduce the observed 
spectral index and its running. We also showed in this Subsection, that when ( is generated 
during inflation and V( is dominated by the tree-level correction, the available parameter 
region is given by Eq. fl3.44p . The other two possibilities are discussed right below. 



C not generated during inflation and V( dominated by the one-loop correction 

According to Eqs. fl4.3p and fl3.33p we have in this case 

= |exp[4iV(r/^-r^.)](^y (^^MfcL), (4.27) 

which reduces to 

{^) » Q'n|exp[4iV(r/^-r/<,)]ln(fcL), (4.28) 
where must be replaced by the observed value V^'^ = (4.957 ± 0.094) x 10~^ |35] . 



( not generated during inflation and P^ dominated by the tree-level term 

The equation 13.441 tells us that in this case the constraint to satisfy is 




(4.29) 



Non-Gaussianity: tnl 



In this section we will calculate the level of non-gaussianity represented in the parameter 
tnl- Since the contributions to the trispectrum T^, calculated in Eqs. f l4.6p and (14. 7p . 
and coming from Figs. IA.5b and IA.5b respectively, present a wavevector dependence as 
that of the first line in Eq. ( I2.32p . we conclude that for the specific subclass of inflationary 
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models we are considering, the non-gaussianity in is parameterized in terms oiTNL- This 
automatically leads to qnl tnl- In view of this, tnl is given in this case as [31] : 



-Tnl 



C 



87r6 



fcifc||k3+k4|3 + 23 permutations 



(4.30) 



The intermediate 0^ T-region 

The level of non-gaussianity t^l according to Eqs. f l4.2p . f l4.7p and fl4.30p . is given by 



1 



l~loop 



+ 23 permutations (Pf-)^ 



||exp[8iV(|r/.|-|r/J)] f ^ V In(fcL) 



.2vr^ 

= ^exp[8iV(|ry.|-|r^,|)](^y^ln(A;L) 

= 2r/^exp[8iV(|r/.|-|r/^|)]Pcln(A;L) 
^ Knl ^ 4.91 X 10^^|r/^|^exp[4001n(5.657 x 10^V^i/2)(|^^| _ 0.020)], (4.31) 

where in the last line we have used the expressions in Eqs. f l3.44p and (13.571) . 

Now, by implementing the spectral tilt constraint in Eq. (13.521) in the spectrum normal- 
isation constraint in Eq. (I3.44p and the amount of inflation constraint in Eq. (I3.58p . we 
conclude that the tensor to scalar ratio r is bounded from below: r > 2.680 x 10~^. 

In the r vs |?7cr| plot in figure l^?T| we show lines of constant tnl corresponding to the values 
TjsfL = 20, 560, 2 X 10^. We also show the high (in white) and intermediate (shaded) 0^ 
T-regions in agreement with the constraint in Eq. (14.231) : 

8 ri4, \fnpj 2 r]^ 

4.070x10^ < |r/^|^exp[4001n(5.657xl0~V-^/2)(|r/<,|-0.020)] < 1.656xl0^^ (4.32) 

As is evident from the plot, the observationally expected 2a range of values for WMAP, 
\tnl\ ^ 2x10^ |112] . PLANCK, \tnl\ ^ 560 |112] . and even the 21 cm background anisotropics, 
ItnlI > 20 and for positive t^vl, are completely inside the intermediate 0* T-region as 
required. Higher values for Tnl, up to t^vl = 1-7 x 10^ are consistent within our framework 
for the intermediate 0^ T-region. 



In subsection 13. 6[ we studied /nl for the case when ( is generated during inflation, 5^ is 
dominated by the one-loop correction, and is dominated by the tree-level term. Fig. 
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Figure 4.1: Contours of tnl in the r vs \ricr\ plot. The intermediate (high) (p^ T- region 
corresponds to the shaded (white) region. The observationally expected 2a range of values, 
for WMAP, PLANCK, and even the 21 cm background anisotropies, and for positive tnl, 
tnl > 20 are completely inside the intermediate 0* T- region. Notice that the boundary line 
between the high and the intermediate 0^ T-regions matches almost exactly the tnl = 0.04 
line. 

[Sil shows the results found. The WMAP [159] (and also PLANCK [I20]) observationally 
allowed 2a range of values for negative /atl, — 9 < /atl, is completely inside the intermediate 
(pi, regioij^. Fig. 14.21 shows both Figs. 13.21 and 14.11 in the same plot. Incidentally, for 
the available parameter window, lines for constant tjvl almost exactly matches lines for 
constant f^i- Thus, it is possible to see that, according the observational status presented 
in the Section l?!4| non-gaussianity is more likely to be detected through the trispectrum than 
through the bispectrum, for the inflationary model studied in this chapter with concave 
downward potential, and from the WMAP, PLANCK, and even the 21 cm background 
anisotropies observations. Fig. 14. 21 also shows some consistency relations between the values 
of Jnl and tnl that will be useful at testing the inflationary model considered with concave 
downward potential against observations. For instance, if WMAP detected non-gaussianity 
through the trispectrum with tnl > 8 x 10^ at the 2a level, the slow-roll inflationary model 




*The intermediate 0* T- region (where T(; is dominated by the one-loop correction and is domi- 
nated by the tree-level term) encloses the intermediate c/j^, region (where is dominated by the one-loop 
correction and P(^ is dominated by the tree- level term). 
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Figure 4.2: Contours of both /atl and tatl in the r vs |?7o-| plot. The intermediate (high) 
0* region corresponds to the shaded (white) region. Lines for constant ttvl almost exactly 
matches lines for constant /atl- According to this figure, and to the observational status, 
non-gaussianity is more likely to be detected through the trispectrum than through the bispec- 
trum, for the inflationary model studied in this chapter with concave downward potential, 
and from the WMAP, PLANCK, and even the 21 cm background anisotropics observations. 
These lines also show some consistency relations between the values of f^i and tnl that 
will be useful at testing the inflationary model considered with concave downward potential 
against observations. 

with concave downward potential considered in this chapter would be ruled out since the 
predicted /tvl would be outside the current observational interval. 

Similarly to the f^i case, it is interesting to see a slow-roll inflationary model with canonical 
kinetic terms where large, and observable, values for tnl may be obtained (in contrast to 
the expected tnl ^ 0{r) from the tree- level calculation |188[ 1192] ). So we conclude that 
if Tc_ is dominated by the one- loop correction but Pc_ is dominated by the tree-level term, 
sizeable non-gaussianity is generated even if ( is generated during inflation. 
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The high 0^ T-region 

According to the expressions in Eqs. ([235]), flH32]) . (K^ . (K^ and (OU]) the value of Tnl 
is in this case 

-j^ r-ptree 

-JNL = p ^ ^ ^ = -r/J = 2 X 10-^ , (4.33) 



87r6 



^i^^^ + 23 permutations (P< 



•>tree\3 



2 



in agreement with the general expectations of Ref. [40] for slow-roll inflationary models 
with canonical kinetic terms where only the tree-level contributions are considered and the 
field perturbations are assumed to be gaussian. This result is no observational interest 
because the generated non-gaussianity is too small to be observable. 



( not generated during inflation 



We will assume in this Section that the fields driving inflation have nothing to do with 
the generation of (; nevertheless, they will generate the primordial non-gaussianity (see for 
instance Refs. [211 M M [M M CDHl [ffl 1201] )• To this end, the post- inflationary 
evolution, particularly the generation of (, will be assumed not to generate significative 
levels of non-gaussianity in comparison with those generated during inflation. 



The low 0^ T-region 

It is possible, in principle, that is dominated by the one-loop correction as long as ( is 
not generated during inflation. Thus, the observed spectral index constraint is no longer 
required and, therefore, the low 0^ T-region is in principle viable. 

Combining the conditions in Eqs. fl4.20p and (14.281) with the expression for the number of 
e-folds in Eq. (I3.57p . we get: 

2 

Tfl 

1 < —V^exp[N\r]^\{4 - 2/n)] , (4.34) 

and 

— ) Vcexp{4N\r^„\), (4.35) 
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where we have defined the parameter n as the ratio between the two rj parameters: n = 
rja/r]^. These two expressions lead to 

rn < 1.6 X 10-^\ri^\ exp{-2N\ri^\/n) , (4.36) 

as a necessary but not sufficient condition to satisfy both Eqs. (14.341) and (14.351) . However, 
by introducing such a condition in Eq. (14.341) . we see that the latter translate into the 
following constraint: 

1 < 10-''|r/^pP^exp[4iV|r/.|(l - 1/n)] . (4.37) 

The previous expression is impossible to satisfy because the highest value the right hand 
side may take is for n — )■ oo and, of course, rj^j = 0.1 and = 62. Such a value, 1.45 x 10~^, 
is much less than one. We conclude that this case is of no interest because it is imposible 
to satisfy the normalisation spectrum condition in Eq. (I4.28p . 



The intermediate 0^, T-region 



The level of non-gaussianity t^tl in this case is given by 

^ rpl — loop 



-TNL 



87r6 



km\^^,+^.4' + 23 permutations 



^"'^ exp[8N{K\ - \vM ] VfHkL) 



^^'exp[8iV(|r/.| - |r/,|)] (^^^ VcHkL) 



exp(8iV|ry.|) - P^ln(A;L) 



^ 2.60 X 10''' (nr)^ (4.38) 

where in the last line we have introduced again the ratio n defined in previous subsubsection, 
and chosen for simplicity |?7o-| =0.1 and = 62 so that the non-gaussianity is maximized. 

From Eqs. (I4.23p . (I4.29P and (I3.57p . and those coming from the ( series convergence 
constraints in Eqs. (I3.76P and (I3.77p . with |?7o-| = 0.1 and = 62, lead to the following 
conditions that reduce the available parameter space: 

• The perturbative regime constraint |x| ^ 1: 

r < 6.51 X 10^ nexp 



12.4 



n 



(4.39) 
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The perturbative regime constraint \y\ <^ 1: 

2.68 X 10-1 



< 



n- 



The dominated by the tree-level term constraint: 



< 



1.10 X 10" 



n 



exp 



12.4 



n 



The Tr dominated by the one-loop correction constraint: 



> 



4.68 X 10" 



-12 



■ exp 



The spectrum normalisation constraint: 



< 



1.6 X 10" 



n 



exp 



37.2 



n 



12.4 



n 



(4.40) 



(4.41) 



(4.42) 



(4.43) 



Analysing these expressions, we conclude that the constraint in the first item is automati- 
cally satisfied once the constraint in the fifth item is satisfied. Moreover, from the constraint 
in the fifth item, we see that the highest possible value r may take is 4.75 x 10~^. And 
finally, to make the constraint in the fourth item consistent with the constraints in the sec- 
ond, third, and fifth items, the lower bound n > 2.58 is required. The resultant available 
parameter window, together with the lines for constant values of tnl, tnl = 1,5,10,15, 
is presented in Fig. lOl for 2.58 < n < 200. Fig. 144] shows the range 200 < n < 2000 
with the lines t^l = 1;5, while Fig. 14.51 shows the range 2000 < n < 3000 also with the 
lines ttvl = 1, 5. As the figures reveal, when is dominated by the one-loop correction and 
is dominated by the tree-level term, large values for t^l are obtained although not so 
large as in the case where ( is generated during inflation. Indeed, an upper bound on t^l, 
according to Eqs. f l4.38p and (14.43 p . is 34.078 when n — )■ oo. 

We conclude that, even if ( is not generated during infiation, we may find observable values 
for tnl. However, such observable values could only be observed by the 21 cm background 
anisotropies at the la level according to the observational status presented in Section 12.41 
We also conclude that, for non- gaussianity to be observable, primordial gravitational waves 
must be undetectable. 
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Figure 4.3: Contours of ttvl in the r vs n plot, for 2.58 < n < 200, when C, is not 
generated during inflation. The allowed parameter space corresponds to the white region. 
The constraint in Eq. f l4.42p almost matches (visually) the horizontal axis. The largest 
possible value tnl may take in this range is 15. 



The high 0^ T-region 



This case is of no interest because the generated non-gaussianity is too small to be observ- 
able: 



1 



n^tree 



87r6 



kiki\i+^.,\^ + 23 permutations 



rptree 



^'ptree^ ' 



_ klki\M4' + 2^ permutations] (Pf-) 
according to the expressions in Eqs. (14.21) and (14. 6p . 



( 'T)tree\ 
2\' C ) 



n 



(4.44) 
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Figure 4.4: Contours of Tnl in the r vs n plot, for 200 < n < 2000, when ( is not 
generated during inflation. The allowed parameter space corresponds to the white region. 
The constraint in Eq. f l4.42p matches (visually) the horizontal axis. The largest possible 
value TjsfL may take in this range is a bit higher than 5. 

1.7.2; fNL 

The low 0^ region 

From Eqs. f l3.34p and (14.201) we see that the low 0^ region and the low (p^ T-region are 
exactly the same, being only constrained by the fact that is dominated by the one-loop 
correction. Thus, the obtained conclusions in Subsubsection UTTT] equally apply. Therefore, 
this case is of no interest because it is impossible to satisfy the normalisation spectrum 
condition in Eq. (I4.28P for the low 0^ region. 
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Figure 4.5: Contours of tnl in the r vs n plot, for 2000 < n < 3000, when ( is not 
generated during inflation. The allowed parameter space corresponds to the white region. 
The constraint in Eq. fl4.42p matches (visually) the horizontal axis. The largest possible 
value tnl may take in this range is a bit higher than 5. 



The intermediate 0^ region 



The level of non-gaussianity is in this case given by 

= ^exp|6JV(|,,,| - (I)Vcln(AL) 

^ -8.59 X 10^(nr)^ (4.45) 

where in the last line we have chosen again for simplicity |?7o-| = 0.1 and = 62 so that 
the non-gaussianity is maximized. 

Since the spectrum normalisation constraint in Eq. f l4.43p equally applies to this case, we 
conclude from it and from Eq. (I4.45P that an upper bound on IJnl] is 2.93 x 10~^ when 
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n ^ oo. fNL is, of course, unobservable. We conclude that when ( is not generated during 
inflation, but the primordial non- gaussianity is, it is impossible to detect non-gaussianity 
through the bispectrum. However, in view of Subsubsection I4.7.H it is possible to detect it 
through the trispectrum. 

The high 0^ region 

This case is of no interest because the generated non-gaussianity is too smaU to be observ- 
able: 

according to the expressions in Eqs. fl3.27p and fl3.29p . 



SECTION 4.8 



Conclusions 



In this chapter we extended the analysis given in the previous one, but this time we cal- 
culated the trispectrum of the primordial curvature perturbation (, generated during a 
slow-roll inflationary epoch and considering a two- field quadratic model of inflation with 
canonical kinetic terms. In order to obtain a large level of non-gaussianity, we consider 
loop contributions as well as tree level terms, and show that it is possible to attain very 
high, including observable, values for the level of non-gaussianity tnl if is dominated 
by the one-loop contribution and is dominated by tree level term. The statement pre- 
sented in Ref. [38] about the suppression of the loop corrections against the tree-level terms 
when considering classicality was analyzed in this chapter and argued to be too strongly 
stated leading to non-general conclusions. The probability that a typical observer sees a 
non-gaussian distribution in the model considered in this thesis was investigated and found 
to be non- negligible. 



Chapter 



NON-GAUSSIANITY FROM VECTOR FIELD 

PERTURBATIONS 



I SECTION 5.1 

Introduction 



The anisotropies in the temperature of the cosmic microwave background (CMB) radiation, 
which have strong connections with the origin of the large-scale structure in the observable 
Universe, is one of hottest topics in modern cosmology. The properties of the CMB tem- 
perature anisotropies are described in terms of the spectral functions, like the spectrum, 
bispectrum, trispectrum, etc., of the primordial curvature perturbation ( [13]. In most of 
the cosmological models the n-point correlators of ( are supposed to be translationally and 
rotationally invariants. However, violations of such invariances entail modifications of the 
usual definitions for the spectral functions in terms of the statistical descriptors [H [121 SI] • 
These violations may be consequences either of the presence of vector field perturbations 

(niiiniiiniEniEiiiHaEsiiaiEsiiHizHii^ mm, spinor 

field perturbations [301 1194] , or p-form perturbations [701 [TH lllll I115[ 1116] , contributing 
significantly to C, of anisotropic expansion [HI |30l SH EH EH HOH [HSl MSI MB l2T8] 
or of an inhomogeneous background [121 El ESj- Violation of the statistical isotropy 
(i.e. violation of the rotational invariance in the n-point correlators of C) seems to be 
present in the data [TH EDI EH 1178] and, although its statistical significance is still low, 
the continuous presence of anomalies in every CMB data analysis (see for instance Refs. 
[32[58l[62l|63lE5lE6l|9ll|92l [1221 11231 [I62l MM 1205]) suggests the evidence might be 
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decisive in the forthcoming years. Since the statistical anisotropy is observationally low, it 
entails a big problem when vector fields are present during infiation, because they generically 
lead to a high amount of statistical anisotropy, higher than that coming from observations 
[52| [75| 1102] . To solve this problem, people use different mechanisms in order to make 
those models consistent with observation, for example using a triad of orthogonal vectors 
[T^ [25], a large number of identical randomly oriented vectors fields [73], or assuming that 
the contribution of vector fields to the total energy density is negligible [521 1102] . 

The amount of statistical anisotropy is quantified throught the parameter (7^, usually called 
the level of statistical anisotropy in the spectrum. Eq. f l2.34p gives us the primordial power 
spectrum that takes into account the leading effects of violations of statistical isotropy by 
the presence of some vector field in the infiationary era. As we could see in Section 12.41 
the gc^ parameter has observational bounds and works, together with the non-gaussianity 
parameters /nl, tnl, S'nl, etc., as statistical descriptors for (. Therefore, it could be a 
crucial tool to discriminate between some of the more usual cosmological models. 

Recent works point out the possibility that a vector field causes part of the primordial cur- 
vature perturbation and show that the particular presence of vector fields in the infiationary 
dynamics may generate sizeable levels of non-gaussianity described by /nl [231 1103[ 1210] 
and tnl [2011209] . In such works the authors included both vector and scalar field perturba- 
tions, and asummed that the contributions to the spectrum from vector field perturbations 
were smaller than those coming from scalar fields and in an opposite way for bispectrum 
and trispectrum. 

In this chapter we use the SN formalism to calculate the tree-level and one-loop contri- 
butions to the bispectrum and trispectrum of (, including vector and scalar field 
perturbations. We then calculate the order of magnitude of the levels of non-gaussianity 
in and including the one-loop contributions and write down formulas that relate the 
order of magnitude of the levels of non-gaussianity /nl and tnl with the amount of statis- 
tical anisotropy in the spectrum gf^. Finally, comparison with the expected observational 
bound from WMAP is done. 



■.SECTION 5.2 

Statistical descriptors from vector field perturbations 



As we saw in Chapter [21 the 5N formalism [521 [1391 ESI [Ml [2021 [m] is a powerful tool 
to calculate the primordial curvature perturbation and all its statistical descriptor to any 
desired order. In the simplest case where ( is generated by one scalar field and one vector 
field and assuming that the anisotropy in the expansion of the Universe is negligible, it can 
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be calculated up to quadratic terms by means of the following truncated expansion [52] : 

C(x) = (5iV(0(x),A,(x),t) = N^6<f)+N\6A + ^N^^{6<j)y+N;j(j)6A+^NjJA6A,, (5.1) 
where 

dN ■ dN d'^N ■• d^N ■ d'^N 

^''^W ^^^^dA' ^"""^W ^^^^^^dlm,' ^^^^Mla^' ^^-^^ 

being the scalar field and A the vector field, with i denoting the spatial indices running 
from 1 to 3. In the Section we defined the power spectrum P^, the bispectrum 5^ and 
trispectrum for the primordial curvature perturbation, through the Fourier modes of ( 
as: 

(C(k)C(k')) ^ (2^)35(k + kOPc(k) ^ {27Tf6{k + k')—rc{k), (5.3) 
(C(k)C(kOC(k")> = (2^)3j(k + k' + k")i?c(k,k',k") 

= (2^)35(k + k' + k")p^fic(k, k', k") . (5.4) 

(C(ki)C(k2)C(k3)C(k4)) = (27r)4<5(ki+k2 + k3 + k4)rc(ki,k2,k3,k4) 

Using Eq. f l5.ip and the definitions given in Eqs. fl5.3p and f l5.4p . it was found in Ref. 
that the tree-level contribution to the spectrum is of the form shown in Eq. fl2.36p . that is 

P^{k) = P-°{k)^l + g^{d.kr) . (5.6) 

In addition, an analogous form for the contribution to /nl was given in Ref. |103j . showing 
that both, and /nl have anisotropic contributions coming from the vector field pertur- 
bation. The one- loop correction to the spectrum was also given in Ref. [52], however they 
kept it in an integral form. In this chapter we give the tree-level and one-loop contributions 
to the bispectrum and to the trispectrum. We also estimate the integrals coming from loop 
corections in order to get an order of magnitude for /nl and Tnl- 

Using the Fourier fmodes for Eqs. Eq. (15. ip . (15. 3p . (15. 4p and (15. 5 p and considering contri- 
butions up to one- loop order, the expressions for P,^, and 7^ , ar^: 

Pf-(k) = N'^Vs^ik) + N\N^;,%,{k) 

= N^Ps^ik) + NlP+{k) + (Na • kfP+{k) (riong - 1) , (5.7) 



(27r)35(ki + k2 + ka + k^) ,^^J,,_ \ ,3 rc(ki, k2, kg, k4) .(5.5) 



pi-i°°P(k) = j 



47r|k + ppp^ 
1 



+ 2^AA<A'7Ifc(k + p)7^7(p 



\Nl^Ps^{\k + v\)Ps<t>{p) + %aKaP5^{\^ + Pl)7;,(p) 

(5.8) 



^ These expresions can be calculated using the diagrammatic tool that we will present in a forthcoming 
paper pOSj . 
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«(k, k', k") = NlN^4Vs4k)Vs4k') + c. p.] + A^AiViiVTl 7I„(k)rfc„(k') + c. p. 



^-'°°P^k,k',k" 



c 



'Pscpik)Tij{k') + 5 perm. 



(5.9) 



47rp3|k + pP|k' - p|3 

kl ATmn I a pK K 



47rp^|k + pp|k' — pp 
d^pk^k'^ 



V5<pip)Vs^{\k + p\)Vs4\k' -p\) 

Til{p)Tkn{k + p)Tjm{k" - p) 

rs4p)rs4\^" + p\)Tij{k'-p) 



+ Vs<,,ip)rs4\k' - p\)Tijik" + p) + rs4\k' - p\)rs4\^" + p|)7I,(p: 

ATi AtJ ATkl f d^pk^k''^ 

+ N^aN^^Naa J 4^p3|k. + p|3|k._p|3 



Vscp{p)Tik(k' -p)Tii{k" + p) 



n^lk" + p|)7;fc(p)7^7(k' - p) + n^dk' - p|)7;fc(p)7^7(k" + p) , (5.10) 



7^'-^'^(ki,k2,k3,k4) = NlNl^[Pscpik2)Vs4k4)Vs4\ki + k2\) + npeim.] 

7ife(k2)75m(k4)7in(ki + k2) + 11 perm. 



+ N^N^N^a^aNTa 



'Ps4,{k-2)Vs^{ki)Tij{ki + k2) + 11 perm. 



+ N\n{n\^N'a^ \Tik{k2)T,i{ki)V5c^{\ki + k2|) + 11 perm. 
+ N^N^^N\n{^ \V5<t>{k2)Tij{ki)V5^{\ki + ksl) + 23 perm. 
+ N^N\N\^NfA \P8ci>{k2)Tik{kA)Tji{ki + ks) + 23 perm.j , (5.11) 



7^y°°P(ki,k2,k3,k4) 



d^p k\kl\k^ + k4| 



NlA^fAN^^^N^^A / i^|k, _ p|3|k3 + p|3|k3 + k4 + p^ 



where 



and 



x7Im(p)75fc(ki - p)7;,p(k3 + p)7io(k3 + k4 + p) , 
7-,(k) = T--{k)V^{k) + ^T^;'\k)V^{k) + T^f ^(k)Pio„g(A:) , 



(5.12) 
(5.13) 



7;7^"(k) = 5,, - , T^f '^(k) = e,,fcfcfc , T'°'^^(k) = kk ■ (5-14) 



Eq. (15. 7p was writen in the form of Eq. (15.61) with d = Na, ^a being a vector with 
magnitude Na = \/ N\N\, and riong = where Piong is the power spectrum for 

the longitudinal component, and P+ and V- are the parity conserving and violating power 
spectra defined by 

V±^]^{Vr±Vl) , (5.15) 

with Vr and Vl denoting the power spectra for the transverse components with right- 
handed and left-handed polarisations [52] . 
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The above expressions can be further separated into different terms: one due to perturba- 
tions in the scalar field, another due to the vector field perturbations, and the other due to 
the mixed terms: 

Pc*=--(k) = Vcf^ik)+V^Ti^), (5.16) 
^i-ioop(j^) = p^i-i-P(A;)+p^^"'°°P(k)+P4^^°°P(k), (5.17) 

Sf-(k,k',k") = ScJ'=^(k,k',k") + i3c7^(k,k',k") + S4T(k,k',k"), (5.18) 
i3j-'°°P(k, k', k") = B^y^^^ik, k', k") + i3ci"'°°P(k, k', k") 

+ B,l'-J^''^{k,k',k"), (5.19) 

7j'-^«(ki, k2, kg, k4) = Tcf'ik,, k2, kg, k4) + rc7^(ki, k^, kg, k4) 

+ r47(ki,k2,k3,k4), (5.20) 

'7^^~'°°''(ki, ka, kg, k4) = 7^J,~^°°P(ki, ks, kg, k4) + 7^^"'°°P(ki, kg, kg, k4) 

+ rci:;^°°P(ki,k2,k3,k4), (5.21) 

Observational analysis tell us that the statistical anisotropy in CMB temperature perturba- 
tion could be observable in a future through current experiments like WMAP or PLANCK. 
Eq. (15.61) combined with recent studies [80] tells us that the level of statistical anisotropics 
gt^ has an upper bound and in the best case (99% confidence level) this is (7^ < 0.383 [75] . 
During our analysis we will adopt an upper bound for (7^: < 0.1. In order to satisfy the 
latter observational constraint over the spectrum, we must be sure that the contributions 
coming from vector fields in Eqs. (15.71) and (15.81) are smaller than those coming from scalar 
fields. That means that the first term in Eq. (15.161) dominates over all the other terms, 
even those coming from one-loop contributions. With the previous conclusion in mind we 
feel free to make assumptions over the other contributions, specially for those coming from 
vector field perturbations. 

^SECTION 5.3" 

Vector field contributions to the statistical descriptors 



As we explain in the previous section, our unique restriction from observation is related to 
the amount of statistical anisotropy present in the spectrum, so we need to be sure that 
the first term in Eq. (I5.16P always dominates. In our study we will assume that the terms 
coming only from the vector field dominate over those coming from the mixed terms and 
from the scalar fields only, except for the case of the tree-level spectrumH. Based on the 



^For an actual realisation of this scenario, we need to show that such constraints are fully satisfied. 
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assumption made, Eqs. f l5.16p - fl5.19p take the form: 




P^^-'°°P(k) 
i3f'='=(k, k',k") 
i3j-'°°P(k, k',k") 

7j^^^(ki,k2,k3,k4) 

^^^°°^(ki,k2,k3,k4) 



rc7^(ki,k2,k3,k4), 

7^\"^°°P(ki,k2,k3,k4) , 



Scr(k,k',k"), 
i3ci-^°°P(k,k',k") 



(5.22) 
(5.23) 
(5.24) 
(5.25) 
(5.26) 
(5.27) 



The above expressions lead us to eight different ways that allow us to study and probably 
get a high level of non- gaussianit}u 

• Vector field spectrum (V(^^) and bispectrum (-B^^) dominated by the tree- level terms 



• Vector field spectrum (P^) and bispectrum (i3^^) dominated by the 1-loop contribu- 



• Vector field spectrum (V(^^) dominated by the tree- level terms and bispectrum (i3^^) 
dominated by the 1-loop contributions. 

• Vector field spectrum (P(^) dominated by the 1-loop contributions and bispectrum 
(i3^^) dominated by the tree-level terms. 

• Vector field spectrum (P^^) and trispectrum (7^^) dominated by the tree- level terms. 

• Vector field spectrum (Vq^) and trispectrum (7^^) dominated by the one- loop contri- 
butions. 

• Vector field spectrum (V^^) dominated by the tree- level terms and trispectrum (7^^) 
dominated by the 1-loop contributions. 

• Vector field spectrum {V(^^) dominated by the 1-loop contributions and trispectrum 
{'Tc.a) dominated by the tree-level terms. 

In order to study these possibilities, we first need to estimate the integrals coming from 



loop contributions. From Eqs. (ES]), CTl . fl512l) . ( 15:231) . (K25\i and (K27\i the integrals 



tions. 



^Our assumption is inspired in the one given in Ref. |31| . 
instead of one scalar and one vector field as in this chapter, 
in Chapter [3](sec also Refs. [H [T7S1 \m\ ). 



In that work the authors use two scalar fields 
A realisation of such a scenario can be found 
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to solve are: 

rc'-'^°n^) = In^IaN'A I ,J!\1%\3 ^^^(^ + P^^dp), (5.28) 

xTil{p)Tkn{^ + p)Tim{^' -p) (5.29) 

x7;^(p)7^fc(ki - p)r„p(k3 + p)7Io(k3 + k4 + p) .. (5.30) 

The above integrals cannot be done analytically, but they can be estimated in the same 
way as that presented in Refs. [31], 11341 1138] . In Appendix [B] we show that the integrals are 
proportional to ln(A;L) where L is the box size. To evaluate them we take the spectrum to 
be scale-invariant, which will be a good approximation if both scalar field (p and vector field 
A are sufficiently light during inffation. The integrals are logarithmically divergent at the 
zeros of the denominator and in each direction, but there is a cutoff at ~ L~^. We found 
that in our case the integrals are also proportional to In(fcL) and that each singularity gives 
equal contributions to the overall result. We find from Eqs. (15.281) and (I5.29p : 

^cY°°'(k) = lNj^NW^V^ + Vion,)S.kr,i{k)\n{kL), (5.31) 
i3j"^°°P(k, k', k") = N'i^N'^^NZ HkL) + Viong)5a [rfc„(k)7-„(k')] (5.32) 

r/^-^°°p(ki, k2, ks, k4) = iv2i<^ivxiiv7^ HkL) {2v+ + p^onJ^^^T-^ki) 

xr.p(k3)7Io(k4 + k3). (5.33) 



^ ^TSdTION 5.4' 



Calculation of the non-gaussianity parameter /] 



NL 



The non-gaussianity parameter in the biespectrum 5^ is defined by [148[ ll2ofi 

/nl = ^ gc(k,k-,k-) _ 

6 [V(;{k)Vi^{k') + eye. perm.] 

Since the isotropic contribution to the curvature perturbation is always dominant compared 
to the anisotropic one, we can write in the above expression only the isotropic part of the 
spectrum P(;^'°(A;): 

f _5 B^{k,k',k") 

/NL g [p^-°(A;)p^i-(A;') + eye. perm.] ' ^ ' ' 



^We employ the WMAP sign convention. 
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Keeping in mind the above expression, we will estimate the possible amount of non- 
gaussianity /nl in the bispectrum B(;, generated by the anisotropic part of the primordial 
curvature perturbation. To do it we take into account the different possibilities mentioned 
in the previous section, where the non-gaussianity is produced solely by vector field pertur- 
bations. 



Vector field spectrum and bispectrum {B(^^) 

dominated by the tree-level terms 

We start our analysis by considering the case studied in Ref. [lU3j . where the authors assume 
that the bispectrum is dominated by vector fields perturbations and that the higher order 
contributions from the vector field are always sub-dominant, i.e N\6Ai 3> N^/j^6Ai6Aj. 
This means that both the spectrum and the bispectrum are dominated by the tree level 
terms, i.e. V(^^^°^ ^ and Bt^^^^^ ^ so that the level of non-gaussianity J'nl 

is given by: 

/NL g ^p^iso^^^p^iso^^,) ^ pg^^ j ~ g [P^^^°(A;)P^'^°(A;') + eye. perm.] ' ^ ' ' 

Since the anisotropic contribution to the curvature perturbation is subdominant, we can 
take ~ 'P(^^^°, so we may write: 

, N\N'^N^X[nMrknik') + eye, perm.] 
YPc_{k)P(_[k') + eye. perm. J 

Assuming that Piong, and V- are all of the same order of magnitude, and that the 
spectrum is scale invariant, we may write the above equation as: 

/..-^^, (5.38) 

where Va = 2P+ + V\ong- Taking as a typical value for the vector field perturbation 
6 A = ^/Va and NaSA > Naa^^"^, the contribution of the vector field to ( is given by 
~ V^Ca ~ ^ as/Va- Thus, we may write an upper bound for /nl: 



-p 3/2 

/nl < (5.39) 



Since the level of statistical anisotropy in the power spectrum is of order gc^ ~ '^QaI'^C^ ^^'^ 
since ^ 5 X ^g-^ [110], Eq. fOSj) yields [IDS]: 

/nl < (l^)'^'- (5-40) 

The above expression gives an upper bound for the level of non-gaussianity /nl in terms of 
the level of statistical anisotropy in the power spectrum when the former is generated 
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by the anisotropic contribution to the curvature perturbation. As we may see, the recent 
observational bounds on /nl: —9 < /nl < HI |119j^ - niay easily be exceeded. 

As an example of this model, we apply the previous results to a specific model, e.g. the 
vector curvaton scenario [49l |50l [51], where the A^-derivatives are |103j : 

= ^r, (5.41) 

Naa = ^r, (5.42) 

where A = |A| is the value of vector field just before the vector curvaton field decays and 
the parameter r is the ratio between the energy density of the vector curvaton field and 
the total energy density of the Universe just before the vector curvaton decay. We begin 
exploring the conditions under which the vector field spectrum and bispectrum are always 
dominated by the tree-level terms. From Eqs. (15. 7p . (15. 9p . (I5.3ip and (I5.32p our constraint 
leads to: 

-PaNI » VlNl^, (5.43) 
VIN^Naa » VlNl^. (5.44) 



. (5.45) 



Thus, it follows that: 

Va <^ ( 

Naa 

We have to remember that in the present case the contribution of the vector field to C is 
given by (a ~ ^/'^Ca ~ NaVPa- Then, the above equation combined with Eqs. (I5.65P and 
(I5.66P leads to: 

r > 2.25 X IQ-^gl^^ . (5.46) 

This is a lower bound on the r parameter we have to consider when building a realistic 
particle physics model of the vector curvaton scenario. 

Finally, from Eq. (15. 38 p . the J'nl parameter in this scenario is given by: 

This is a consistency relation between Jnl, QCi ^"^^ which will help when confronting the 
specific vector curvaton realisation against observation. 



^The bispectrum (trispectrum) in this scenario might be either of the local, equilateral, or orthogonal 
type. We are not interested in this thesis on the shape of the non-gaussianity but on its order of magnitude. 
Being that the case, comparing with the expected bound on the local fyh |119j (t>jl |112| ) makes no sensible 
difference under the assumption that the expected bounds on the equilateral and orthogonal /nl ("nsTL) are 
of the same order of magnitude, as analogously happens in the /nl case for single- field inflation |193j . 
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Vector field spectrum and bispectrum {B(^^) 

dominated by the 1-loop contributions 

Since the bispectrum is dominated by 1-loop contributions and is given by Eq. fl5.32p . we 
may write Eq. fl5.35l) as: 

NIaNXaNTa HkL) {2V+ + Viong)Sa [%n{^)r,m{^') + cyc. perm.] 

- r-n n \T> I ^ • {^■^^) 

[P(;{k) P(^[k' ) + cyc. perm. J 

Assuming again that Piong; and V- are all of the same order of magnitude, and that 
the spectrum is scale invariant, the above equation leads to: 



/nl - . (5.49) 



Since the vector field spectrum is dominated by the 1-loop contribution, (a ~ \/'^c,a ~ 
^ aaPa- Thus, and taking into account that gc_ ~ '^c,aI'^(, ^^'^ Vc^^"^ ^ 5 x 10"^ |119j . we 
find: 

The biggest difference between the result found in Ref. |1U3] . given by Eq. fl5.40l) . and the 
result given by Eq. fl5.50p . is that the latter gives an equality relation between the non- 
gaussianity parameter /nl and the level of statistical anisotropy in the power spectrum gc^. 
Following the recent bounds for /nl: —9 < /atl < 111 |119j . this scenario predicts an upper 
bound for the gc^ parameter: 

gc^ < 0.02. (5.51) 
This bound is stronger than that obtained from direct observations in Ref. [80] . 

Again we apply our result to the vector curvaton scenario. Since we are assuming that the 
vector field spectrum and bispectrum are dominated by 1-loop contributions, we get from 
Eqs. ([52D, dtnD, f lCT]) . and fO^ : 

rA>(^]\ (5.52) 



Naa 

which for the vector curvaton scenario becomes: 

r < 2.25 X 10"\^/^ . (5.53) 

This is an upper bound on the r parameter we have to consider when building a realistic 
particle physics model of the vector curvaton scenario. 
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Vector field spectrum ("P^^) dominated by the tree- level 
terms and bispectrum {B(^^) dominated by the 1-loop 
contributions 

In order to check the viabihty of this case, we start studying the imphcations of the restric- 
tions over the spectrum and the bispectrum, i.e. what happens when we assume that the 
vector field spectrum is dominated by the tree-level terms and the bispectrum is dominated 
by the 1-loop contributions. From Eqs. (15. 7p . ( 15. 9p . (I5.3ip . and (I5.32p it follows that: 

VaNI:>VInI^ VA-^jjf, (5.54) 
VInINaa<^V\NI^ VA^jrf. (5.55) 

^^AA 

As we may see, it is impossible to satisfy simultaneously Eqs. ( I5.54p and (I5.55p . This is 
perhaps related to the fact that we have taken into account only one vector field. Such a 
conclusion may be relaxed if we take into account more than one vector field, as analogously 
happens in the scalar multi-field case |l3l 1175] . 



p.4.4 Vector field spectrum dominated by the 1-loop 

contributions and bispectrum (^B^^) dominated by the 
tree-level terms 



As in the previous case, it is impossible to satisfy the conditions under which the spectrum 
is always dominated by the 1-loop contributions and the bispectrum is always dominated 
by the tree-level terms: 



Nl 

» (5.56) 

^^AA 
Nl 

Va « -f. (5.57) 

^^AA 



Again, the conclusion may be relaxed if we take into account more than one vector field. 
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Calculation of the non-gaussianity parameter tnl 



The non-gaussianity parameter tnl in the trispectrum Tr is defined by [3l 



^ _ 2rc(ki,k2,k3,k4) 

"""^ [Pc(ki)n(k2)Pc(ki + k4) + 23perm.] ■ ^' ^ 

Remember that the isotropic contribution in the Eq. fl5.6p is always dominant compared to 
the anisotropic one so that we may write in the above expression only the isotropic part of 
the spectrum 

^ ^ 27^(ki,k2,k3,k4) 

"""^ [Pc"°(^i)n"°(^2)n"°(|ki + k4|) + 23perm.] ■ ^' ' 

Using the above expression, we will estimate the possible amount of non-gaussianity gener- 
ated by the anisotropic part of the primordial curvature perturbation, taking into account 
different possibilities and assuming that the non-gaussianity is produced solely by vector 
field perturbations. 



^.5.1 Vector field spectrum (V(^) and trispectrum (7(^) 
dominated by the tree-level terms 

In this first case, we assume that the trispectrum is dominated by vector field perturbations 
and that the higher order terms in the 6N expansion in Eq. fl5.ll) involving the vector field 
are sub-dominant against the first-order term: N\SAi ^ N^J^^SAiSAj. The latter implies 
that both the spectrum and the trispectrum are dominated by the tree-level terms, i.e. 
VcT > ^Ca^^""" and > TcV^""". Thus, we have from Eq. (l539|) : 

^ ^ 2rcr(ki,k2,k3,k4) 

"""^ [Pr{ki)Vrih)Vri\^i + ^A\) + 23peTmy ^ 



which, in view of Eqs. (15. lip and f l5.26p . looks like: 



Tik{k2)Tjm(Mi)Tin(ki + ks) + 11 perm. 



2N\N^^N^2^N^2 [7-fc(k2)7-^(k4)7I„(ki + k2) + 11 p 
Vr{h)Vr{h)Vr{\^i + k4|) + 23 perm.] 



(5.61^ 



We will just consider here the order of magnitude of tnl- Therefore, we will ignore the 
specific k dependence of Tij- Instead, we will assume that Piong, and V- are all of 
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the same order of magnitude, which is a good approximation for some specific actions (see 
for instance Ref. |52]), and take advantage of the fact that the spectrum is almost scale 
invariant |119] . Thus, after getting rid of all the k dependences, the order of magnitude of 
tnl looks hke: 

Tnl ^ A A AA 5_g2 

'■^^ iso\3 ' V ^ 

where Va = 2P+ + "Piong- Employing our assumption that NaSA > NaaSA'^, and since 
the root mean squared value for the vector field perturbation 6A is \/Va, the contribution 



of the vector field to ( is given by (a ~ y "^Ca ~ NavRa- An upper bound for tnl is 
therefore given by: 



Tnl < (5.63) 

"1^ r-^ Iq-y lSO\3 \ ' 

Since the order of magnitude of gc^ is V(^^/V(^^^° , under the assumptions made above we get: 

T-NL < 8x10'(|y)'' (5-64) 

where {Vl''''f/^ ~ 5 x lO'^ [119] has been used. Eq. (ESI) gives an upper bound for 
the level of non-gaussianity tnl in terms of the level of statistical anisotropy in the power 
spectrum gf^ when the former is generated by the anisotropic contribution to the curvature 
perturbation. Comparing with the expected observational limit on tnl coming from future 
WMAP data releases, tnl ~ 2 x 10^ |112] . we conclude that in this scenario a large level 
of non-gaussianity in the trispectrum of C, is possible, leaving some room for ruling out 
this scenario if the current expected observational limit is overtaken. 

As an example of this scenario, we apply the previous results to a specific model, e.g. the 
vector curvaton model |l9l |50l [51] , where the iV-derivatives are |103] : 

Na = ^r, (5.65) 

Naa = ^r, (5.66) 

where A = |A| is the value of vector field just before the vector curvaton field decays and 
the parameter r is the ratio between the energy density of the vector curvaton field and the 
total energy density of the Universe just before the vector curvaton decay. 

First, we check if the conditions under which the vector field spectrum and trispectrum are 
always dominated by the tree-level terms are fully satisfied. From Eqs. (15. 7p . (15. lip . (I5.3ip 
and (I5.33P our constraint leads to: 

VaNI » PliVl^, (5.67) 
rlNlNl^ » V'^NU, (5.68) 



which mean that the if the vector field spectrum is dominated by the tree-level terms so is 
the vector field trispectrum. An analogous situation happens when the vector field spectrum 
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is dominated by the one-loop terms: the vector field trispectrum is also dominated by this 
kind of terms. As a result, it is impossible that simultaneously the vector field spectrum is 
dominated by the tree-level (one-loop) terms and the vector field trispectrum is dominated 
by the one-loop (tree-level) terms. Following Eq. f l5.67p . we get: 

(5.69) 

which, in view of (a ~ \/'Pca ~ ^ asTPa and Eqs. (15.651) and (I5.66p . reduces to: 

r > 2.25 X 10-\^/^ (5.70) 

This lower bound on the r parameter has to be considered when building a realistic particle 
physics model of the vector curvaton scenario. 

Second, looking at Eq. (I5.62p . we obtain the level of non-gaussianity tnl for this scenario: 

2x 10-^ (g^\i 
Tnl ^ 5 77T • (5-71) 



H VO.l 

This is a consistency relation between tnl, fl'c ^'^'^ which will help when confronting 
the specific vector curvaton realisation against observation. Indeed, a similar consistency 
relation between /nl and gc^ was derived for this scenario in 15.471 |21U] : 

Thus, in the framework of the vector curvaton scenario, the levels of non-gaussianity /nl 
and Tnl are related to each other via the r parameter in this way: 

^NL ^ ^/nl . (5-73) 

in contrast to the standard result 

^NL = ^/^L , (5.74) 
for the scalar field case (including the scalar curvaton scenario) found in Ref. |40j . 



fe.5.2 Vector field spectrum and trispectrum (7^^) 

dominated by the one-loop contributions 

From Eqs. fl5.33p and (15.590 we get 
""""^ [Vr{ki)Vr{h)Vr{\k, + k,\) + 23 perm.] ' ^ ' ' 
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Assuming again that Piong, and are all of the same order of magnitude, and that 
the spectrum is scale invariant, we end up with: 

-NL-^. (5.76) 



/nl ~ W'(^Y^\ (5.78) 



Performing a similar analysis as done in the previous subsection, but this time taking 
into account that the vector field spectrum is dominated by the one-loop contribution and 
therefore (a ~ a/^Ca ~ ^aaPa, we arrive at: 

-NL ~ ~ 8 X 10« (5.77) 

The above result gives a relation between the non-gaussianity parameter tnl and the level 
of statistical anisotropy in the power spectrum qq. 

Now, we call a similar result that we found for the non-gaussianity parameter /nl in 
Eq. ( I5.50p . that is: 

.0.1 

By combining Eqs. fl5.77p and (15.781) we get: 

TNL ~ 8 X IOVnl , (5-79) 

which gives a consistency relation between the non-gaussianity parameters /nl and tnl for 
this particular scenario. The consistency relations in Eqs. (I5.77p . (I5.78p . and (I5.79P will put 
under test this scenario against future observations. In particular, the consistency relation 
in Eq. (I5.79P differs significantly from those obtained when ( is generated only by scalar 
fields (see e.g. Eq. (jEZlD and Ref. |10]). 

Again when we apply our result to the vector curvaton scenario, we get from Eqs. (15. 7p . 
(1511D . (ICTl) . ( 15331) . (I5J5D and (Km : 

r < 2.25 X 10"^^/^ , (5.80) 

which is an upper bound on the r parameter that must be considered when building a 
realistic particle physics model of the vector curvaton scenario. 



As happens in Sections 15.4.31 and I5.4.4[ the other two possibilities are not viable because it 
is impossible to satisfy simultaneously that the vector field spectrum (P^^) is dominated by 
the tree- level terms and the trispectrum (7(^) is dominated by the one- loop contributions, 
or the vector field spectrum {V(^^) is dominated by the one- loop contributions and the 
trispectrum (7^^) is dominated by the tree-level terms. 
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— I section"^ 
Conclusions 



We have studied in this chapter the order of magnitude of the levels of non-gaussianity 
/nl and tnl in the bispectrum 5^ and in the trispectrum T^, respectively, when statistical 
anisotropy is generated by the presence of one vector field. Particularly, we have shown 
that it is possible to get an upper bound on the order of magnitude of /nl if we assume that 
tree level contributions on and domiante over all other terms, this result given in 
the Eq. fl5.40p . On the other hand if we assume that the 1-loop contributions dominate 
over the tree- level terms in both the vector field spectrum (P^^) and the bispectrum (-8^^) a 
high level of non-gaussianity /nl is obtained, /nl is given in this case by Eq. f l5.50p . where 
we may see that there is a consistency relation between /nl and the amount of statistical 
anisotropy in the spectrum g,^. We also have shown that it is possible to get an upper bound 
on the order of magnitude of tnl if we assume that the tree-level contributions dominate 
over all higher order terms in both the vector field spectrum {V(^^) and the trispectrum 
(7^^); this bound is given in Eq. f l5.64p . We also found that it is possible to get a high level 
of non-gaussianity tnl, easyly exceeding the expected observational bound from WMAP, 
if we assume that the one-loop contributions dominate over the tree-level terms in both 
the vector field spectrum (V(^) and the trispectrum (7^^). tnl is given in this case by 
Eq. ( I5.77p . where we may see that there is a consistency relation between the order of 
magnitude of tnl and the amount of statistical anisotropy in the spectrum g^. Two other 
consistency relations are given by Eqs. f l5.78p and fl5.79p . this time relating the order of 
magnitude of the non-gaussianity parameter /nl in the bispectrum with the amount of 
statistical anisotropy (7^ and the order of magnitude of the level of non-gaussianity tnl in 
the trispectrum T^. Such consistency relations let us fix two of the three parameters by 
knowing about the other one, i.e. if the non-gaussianity in the bispectrum (or trispectrum) 
is detected and our scenario is appropriate, the amount of statistical anisotropy in the 
power spectrum and the order of magnitude of the non-gaussianity parameter tnl (or /nl) 
must have specific values, which are given by Eqs. fl5.78p (or fl5.77p ) and fl5.79p . A similar 
conclusion is reached if the statistical anisotropy in the power spectrum is detected before 
the non-gaussianity in the bispectrum or the trispectrum is. 



Chapter 



CONCLUSIONS 



Observational cosmology is in its golden age: current satellite and balloon experiments are 
working extremely well [901 1159] . dramatically improving the quality of data |119j . More- 
over, foreseen experiments [591 1206] will take the field to a state of unprecedent precission 
where theoretical models will be subjected to the most demanding tests. Given such a 
state of affairs, it is essential to study the higher order statistical descriptors for cosmolog- 
ical quantities such as the primordial curvature perturbation (, which give us information 
about the non-gaussianity and about the possible violations of statistical isotropy in their 
corresponding probability distribution functions. 

The slow-roll class of inflationary models with canonical kinetic terms are the most popular 
and studied to date. Inflationary models of the slow-roll variety predict very well the 
spectral index in the spectrum of ( but, if the kinetic terms are canonical, they seem to 
generate unobservable levels of non-gaussianity in the bispectrum and the trispectrum 
of ( making them impossible to test against the astonishing forthcoming data. Where 
does this conclusion come from? The answer relies on careful calculations of the levels of 
non-gaussianity /tv^ and t^l by making use of the SN formalism [22 |I188[ [211(1226] . In this 
framework, ( is given in terms of the perturbation 6N in the amount of expansion from the 
time the cosmologically relevant scales exit the horizon until the time at which one wishes 
to calculate (. 

Due to the functional dependence of the amount of expansion, ( is usually Taylor-expanded 
(see Eq. (13. 3p ) and truncated up to some desired order so that f^L and t^vl are easily 
calculated (see for instance Eq. (13. 7p ). Two key questions arise when noting that it is 
impossible to extract general and useful information from the ( series expansion in Eq. (13. 3 p 
until one chooses a definite inflationary model and calculates explicitly the derivatives. 



82 



First of all, when writing a general expression for f^L or t^l in terms of the N derivatives, 
how do we know that such an expression is correct if the series convergence has not been 
examined? Moreover, if the convergence radius of the ( series is already known, why is each 
term is the ( series supposed to be smaller than the previous one so that cutting the series 
at any desired order is thought to be enough to keep the leading terms? Nobody seems to 
have formulated these questions before and, by following a naive line of thinking, /tvl and 
tnl were calculated for slow-roll inflationary models with canonical kinetic terms without 
checking the ( series convergence and keeping only the presumably leading tree-level terms 
[221 [TM 11881 flOTlimi [226] . 

These two questions have been addressed in this thesis (see Chapt. 3 and 4) by paying 
attention to a particular quadratic small-field slow-roll model of inflation with two com- 
ponents and canonical kinetic terms (see Eq. fl3.17p ). Although the non-diagrammatic 
approach followed in Section 13.71 to find the necessary condition for the convergence of the 
C series in our model might not be applicable to all the cases, we have been able to show 
that not being careful enough when choosing the right available parameter space could 
make the ( series, and therefore the calculation of f^i and t^l from the truncated series 
(e.g. Eq. (13. 7p ). meaningless. We also have been able to show in our model that the 
one- loop terms in the spectrum P^, the bispectrum and trispectrum T,^ of ( could be 
bigger or lower than the corresponding tree-level terms, but are always much bigger than 
the corresponding terms whose order is higher than the one-loop order. If is dominated 
by the one-loop correction but is dominated by the tree-level term, sizeable and observ- 
able values for f^L o,re generated, so they can be tested against present and forthcoming 
observational data, a similar conclusion was reached when the trispectrum is dominated by 
one-loop corrections and the is dominated by the tree- level term. Finally, if both P(j 
and or are dominated by the tree- level terms, /tvl or tnl are slow-roll suppressed (see 
Eqs. 13.631 and I4.33P as was originally predicted in Refs. [2 ^ I21H 1226] . What these results 
teach us is that the issue of the ( series convergence and loop corrections is essential for 
making correct predictions about the statistical descriptors of ( in the framework of the 6N 
formalism, and promising for finding high levels of non-gaussianity that can be compared 
with observations. 

The above disccusion about ( was made assuming that the n-point correlators of C are tran- 
lationally and rotationally invariant. However as we could see in the section [2.4.2[ violations 
of the translational (rotational) invariance (i.e. violations of the statistical homogeneity 
(isotropy)) seem to be present in the data [62l[63l[85l[86l[9ll[92] ([Il|79l|80l|88l UTS]); 
therefore it is pertinent to study theoretical models that include those violations. This 
is the reason why in the chapter [S] we studied the statistical descriptors for ( for models 
with vector field perturbations, which are responsible of violations of statistcal isotropy. 
We studied in that chapter the order of magnitude of the levels of non-gaussianity /nl 
and Tnl in the bispectrum P^ and in the trispectrum T^, when statistical anisotropy is 
generated by the presence of one massive vector field. We have shown that it is possible 
to get an upper bound on the order of magnitude of /nl (see Eq. f l5.40p ) and tnl (see 
Eq. f l5.64p ) if we assume that the tree- level contributions dominate over all higher order 
terms in both the vector field spectrum (P(^^), the bispectrum (-B,^^) and trispectrum (7^^). 
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We also show that it is possible to get high levels of non-gaussianity /nl and tnl, easily 
exceeding the expected observational bounds from WMAP, if we assume that the one-loop 
contributions dominate over the tree- level terms in both the vector field spectrum (P,^^) 
and the bispectrum (-B^^) or in both the vector field spectrum {V(^) and the trispectrum 
0~Ca)- could see that there are a consistency relations between the order of magnitude 
of /nl and the amount of statistical anisotropy in the spectrum g^^ [Eq. fl5.50p ] and between 
the order of magnitude tnl and (7^ [Eq. f l5.77p ]. Two other consistency relations are given 
by Eqs. fl5.78p and fl5.79p . this time relating the order of magnitude of the non-gaussianity 
parameter /nl in the bispectrum 5^ with the amount of statistical anisotropy g,^ and the 
order of magnitude of the level of non-gaussianity tnl in the trispectrum T^. Such consis- 
tency relations let us fix two of the three parameters by knowing about the other one, i.e. 
if the non-gaussianity in the bispectrum (or trispectrum) is detected and our scenario is 
appropriate, the amount of statistical anisotropy in the power spectrum and the order of 
magnitude of the non-gaussianity parameter tnl (or /nl) must have specific values, which 
are given by Eqs. fl5.78p (or f l5.77p ) and fl5.79p . A similar conclusion is reached if the 
statistical anisotropy in the power spectrum is detected before the non-gaussianity in the 
bispectrum or the trispectrum is. 



A 



TREE-LEVEL AND ONE-LOOP DIAGRAMS FOR 
Pc, AND : SCALAR FIELDS 



We show in this appendix the mathematical expressions for the tree-level and one-loop 
Feynman-like diagrams associated with the spectrum P^, the bispectrum an trispectrum 
of C, following the set of rules presented in Ref. [39] • To this end we have taken into account 
the N derivatives for our small- field slow- roll model given in Eqs. fl3.24p . f l3.25p . and f l3.26p . 
After presenting the mathematical expressions, we will estimate the order of magnitude of 
each diagram in order to determine the respective leading terms at tree-level and one-loop 
for both and Bq. 



Figure A.l: Tree-level Feynman-like diagram for P^. The internal dashed line corresponds 
to a two-point correlator of field perturbations. 



Tree-level diagram for Pq 
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I SECTION a!^ 



Tree-level diagram for 



Looking at Fig. \A.1\ we see that P^^'^^ is given by 



Tjtree 



(A.l^ 



Of course, there is only one tree-level diagram for and therefore Eq. (lA.ll) is the associ- 
ated leading tree- level term. 

Our calculation in this appendix goes up to the one-loop diagrams so, in order to have 
complete consistency in the calculation |186j . we should also take into account the one- 
loop correction to the two-point correlator in the field perturbations when calculating the 
diagram in Fig. lA.ll Such a correction has been studied in Refs. |185l 11961 11971 12211 1222] 
where the most general result for single-field slow-roll inflation with Ntotai not very much 
bigger than 62 is |185] 



l~loop 



27r2 



El 

27T 



( 



\27rmp J 



3^ 

-HkL)+f3 
6 



(A.2) 



where L is the infrared cutoff for a minimal box [271 1138] , and /3 is a renormalisation 
scheme-dependent constant that is expected to be negligible on large scales compared to 
ln(/cL) ~ 0{1). The one-loop correction to the field perturbation spectrum in Eq. (lA.2p 
is, therefore, negligible compared to the tree-level contribution P^lp"^ = (27r^ /k^){H^/27r)'^ 
if Hi, ^ mp as usually required. In our model <^ mp is indeed given but, since we 
are dealing with a two-component model, the result in Eq. (IA.2P may not be applicable. 
Anyway, we feel quite confident that the (up to now unknown) extension of Eq. flA.2l) to 
the multiple-field case will yield similar results, so we will keep the expression in Eq. (lA.ip 
as the leading tree- level contribution to P^. 



One-loop diagrams for P,^ 
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(b) 



Figure A. 2: One-loop Feynman-like diagrams for P^. (a). The two internal dashed lines 
correspond to two-point correlators of field perturbations, (b). The internal dashed lines 
correspond to a three-point correlator of field perturbations. 



One-loop diagrams for 



Looking at Figs. IA.2b and IA.2b . we see that P} is given by two contributions P} 



andP^'-'°°^': 



l—loop a 



c 



P, 



l~loop b 



c 



(A.3) 



f d'^q 

J {2it)^ ^^'^ ^'^ S^{k,q, |k + q|) + 
f d^Q 

+N^N„a- J ^^^^50 5<7 |k + q|) 

J (27r)3 



+ 



+ 



.3^3 



exp [2N{r]^ - Va)] 



J (27r)3 



perm. I2a. 



27ry 2^/2mp A;3g3|k_^q|3 



(A.4) 



where the ln(/cL) ~ C^(l) factor comes from the evaluation of the momentum integrals in 
a minimal box [271 1114[ 1138] . the A^3(/ci, ^2, Ajs) function is defined by |189] 



M3{ki,k2,k3 



3j,3 



^3^ 



2'^3 



{k2 - ^3 



(A.5) 



with kt = ki + k2 + k-s, and the subindex perm. I2a. means a permutation over the last two 
arguments in A^3. 



A quick glance reveals that the first term in Eq. ( lA.SP is subleading with respect to the 
second one because |?7o-| > |?70| and exp[4A^(?7<^ ~ Va)] ^ 1- The same is true for Eq. ( ]A.4p 



One-loop diagrams for P,^ 
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where exp[2A^(?7<^ — rja)] ^ 1. Now, by comparing the orders of magnitude of the leading 
terms in Eqs. flA.SP and (IA.4p . we conclude that: 



ryl — loop b . 1/2 

Vamp roAr/^ „ M J- 

e 



exp [2iV(r/^ - r/,)] ^ > 1 , (A.6) 



where mp ^ (p^ and <^ 1. Thus, the one-loop leading term for in our model is given 

by 



Having presented the leading tree-level and one-loop contributions to in Eqs. ( lA.ip and 
(]A.7p . a consistency issue to think about is the dependence of the expression in Eq. ( ]A.2p 
on the infrared cutoff L. This quantity is in principle an artefact of the series expansion, 
and the final series result should in principle be independent on the chosen value for L (see 
for instance Ref. |170] ). In fact, by assuming that this is the case, Refs. [211 EOl 1138] have 
shown that there is a running on the N derivatives with respect to L so that changes in 
the In(fcL) factors are compensated by the running of the N derivatives. This is similar to 
what happens in Quantum Field Theory where physical results independent on the energy 
scale must be independent of the chosen value for the renormalisation scale Q. Changing 
Q only modifies the relative weight of the tree-level and loop contributions, usually making 
the tree- level terms dominate over the loop corrections if Q is chosen around the relevant 
energy scale of the process studied. Nevertheless, we see that the ln{kL) term in Eq. ( ]A.2|) 
does not compensate for the ln(A;L) term in Eq. (lA.7p . which is a real concern as we could 
expect since ( and its spectral functions are a set of observables. The solution to this 
paradox relies on the fact that the observed ( depends on L as the stochastic properties 
of the distributions depend on the size of the available region in which we are actually 
able to perform observations. In this regard ( is analogous to for instance the fine structure 
constant in Quantum Field Theory which, being an observable, depends on the energy scale 
for which experiments are done and, therefore, on Q. Likewise, ( and its spectral functions, 
though being observables, depend on the size of the regions where observations are done 
and, therefore, on L. Having this in mind it is essential to work in a minimal box [21], i.e. 
with L a bit bigger than H^^ (with the subscript meaning today), so that In(fcL) ~ 0{1) 
as has been done throughout this thesis. 



Tree-level diagrams for Bq 
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Figure A. 3: Tree-level Feynman-like diagrams for Bq. (a). The two internal dashed lines 
correspond to two-point correlators of field perturbations, (b). The internal dashed lines 
correspond to a three-point correlator of field perturbations. 
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Tree-level diagrams for Bt 



Looking at Figs. IA.3b and lA.Sb . we see that B^^'^^ is given by two contributions B^^'^^ " and 

c 



j^tree b. 



^tree a ^ ^2^^^ [Ps^{h) Ps^ih) + 2 pemiutatious] 



j^tree b ^ N^Bs^ 54, 5<p{kl, ^2, ^3 



Now, from comparing the order of magnitude of the expressions in Eqs. (1A.8I) and (1A.9I) . 
we conclude that: 



jDtree b . ,1/2 / _ „ 



^ (A.0) 

which in fact is usual as demonstrated in Refs. |147[ 1211] . Thus, the tree- level leading term 
for B(^ in our model is given by: 
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As was done for in Subsection lA.ll the one-loop correction to the spectrum of the field 
perturbations must be taken into account for the sake of consistency when calculating the 
contribution associated to the diagram in Fig. IA.3b . The discussion about the relevance of 
this quantum one-loop correction is actually the same as in Subsection lA.ll and, therefore, 
we may conclude with some confidence that the expression in Eq. flA.Sp is reliable. As 
regards the diagram in Fig. IA.3b . it is necessary to include the one-loop correction the three- 
point correlator of the field perturbations in Eq. (lA.Op . which in fact nobody has calculated 
yet even for the single-field case. Nevertheless we might conjecture that, analogously to 
that for the case, such a correction is negligible compared to the tree-level contribution 
to and, therefore, the expression in Eq. (1A.9P will also be reliable. 
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(a) 



(b) 



(d) 



(e) 



(f) 



Figure A. 4: One-loop Feynman-like diagrams for i?^. (a) and (b). The three internal 
dashed lines correspond to two-point correlators of field perturbations, (c), (d), and (e). 
The internal dashed lines correspond to a two-point and a three-point correlator of field 
perturbations, (f). The internal dashed lines correspond to a four-point correlator of field 
perturbations. 
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I SECTION A!^ 



One-loop diagrams for 



Looking at Figs. IA.4k . IA.4b . IA.4b . IA.4l i. IA.4b . and IA.4f . we see that 5^'"''^ is given by six 
contributions ^i~ioopb^ ^i-ioopc^ ^i-ioopd^ ^i^ioope^ 



r)l— loop a 



ryl—loop b 



Tyl—loop c 



[Nt^ + NL] I 7^n<^(gm<^(|ki + q|)n^(|k3-q|) 



(2vr)^ 



1 Tl"^ 



(A.12) 
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1 
2 



(27r)3 

2 27^3 



-P59i(g)-P5</,(|k3 - q\)Psct>{k2) + 5 permutations 



3^6 ^6^6 - ^-)] 



IGTr"^ ( ) ln(A;L) 



(A.13) 



ks + q|, A;3)P5^(|ki - q|) + 2 permutations] + 



(2^ 



[55a 5a 5</.(g, Iks + q\,k3)Ps^{\k:L - q|) + 2 permutations] 



(fq 



in 



E 



perm 



271 J 2V2mp q^\k3 + q\^ki |ki - q| 



+ 



-2 permutations 



5^5 



exp [4:N{r]^ - Va)] x 



X 



d'q g 6 (H.y ■M3(fc3,g,|k3 + q|) 1 

(2^)' pe.tl2aA2vry' 2v^mp A:3V|k3 + qP |ki - q|3 



+2 permutations 



(A.14) 
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T^l—loop d 



B 



l—loop e 



(2^ 



<5</, 5<pih, q, Iks - q\)P5cf>{k2) + 5 permutations] + 



(27r) 



(^3, g, Iks - (\\)P5^{k2) + 5 permutations] 



.3^5 



(2^ 



57r 



^^2ti) 2V2mp A:|g3|k3-q|3 /i;3 



perm 



+5 permutations 
-^^exp[2iV(r/^-r/„)] x 
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^ /'^V e^' 7^3(^:3, g, Iks -q|) 1 



(2-)^^" ,e.^.,..V2.y 2V2mp fc3V|ks-qP 



-5 permutations . 



(A. 15) 



+-N'iN 



d^q 
d^q 



(27r)3 



[^50 50 5(^{ki, q, |ki + q|)P5<^(/c2) + 5 permutations] + 
[Bsip Sa Scrih, q, |ki + q\)Ps^{k2) + 5 permutations] 



l^*y •^3(fci,g,|ki + q|) 1 ^ 



perm 



271 J 2V2mp /i;?g3|ki + q|3 



+5 permutations — 



5^5 



exp [2N{r]^ - Va)] x 



-1/2 



H.y V 7W3(A;i,g,|ki + q|) 1 



d^q 6 ( 

(27r)3 V27r; 2v^mp A:iV|ki + q|3 kl 

perm. I2a. ^ ' " ^ 



+5 permutations . 



(A.16) 
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Tjl-loop f 



3^4 



(27r)3 
d^q 



[Tscj, Sep 5<p 5<p{^i, q, ka - q, k2) + 2 permutations] + 
[T54, 5a 5a 50 (ki, q, ks - q, k2) + 2 permutations] 



(27r) 



+2 permutations 



Stt' 



+ 



perm 



A^4(ki,q,k3 - q, ka) 1 
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exp [2N{r}^ - r]^)] x 
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^ /^^y -M4(ki,k2,q,k3-q) 1 



perm. /2a. Z2a. 

+2 permutations , 



27r / /c^/cfg^lks — q| 



(A.17) 



where the subindex perm. /2a. I2a. means a permutation over the first two arguments and 
simultaneously over the last two arguments in A^4(ki, ka, ks, k4) defined by |19U] 



A^4(ki,k2,k3,k4) 



Ta\2 ■ Z34 01 rv 

"2 V 2k2 • Z34 + -0'i2Cr34 

^34 4 



1 fci 
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2 ^34 
34 



ki ■ ks 



Vri24 + 



with kjj = kj + kj, = /ci + ^2 + /cs + ^4, and 
Cjj = kj ■ kj + kj , 



Imn 



1 + 
= 1 



1 _^ "^kikj 



kt 

kl + fcm + kn 



2{kikm + /l^i/^n ~l" kmkn) ^ Qkik^nkn 



kf 



(A.18) 

(A.19) 
(A.20) 

(A.21) 
(A.22) 



Following the same kind of analysis as we carried out for the one-loop diagrams of and 
the tree-level terms for we conclude the following: 



j^l— loop a 



B, 



l—loop b 

c 



B 



l—loop a 

_c 

Tjl—loop c 



^exp[4iV(r,,-r,.)]4vr^ (0) (f)^ 
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^ exp [4Ar(^^ _ ^^)] 4vr4 !r 



Va mp 
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1 
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(A.23) 



(A.24) 
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rsj 

■Til — loop d ^ ^1/2 ^„ 3 

#f exp [2N{^, - ^.)] (f )^ ^ 
^"j ^exp[4iV(r/^-r/.)]4i»l, (A.25) 



e 



- — exp [47V(r/^ - r/^)] ^ » 1 , (A.26) 



Ej-'-^"^ ^llf exp [6iV(r/<^ - r/.)] (g^) (t)'47r4 



Thus, the one-loop leading term for i?^ in our model is given by: 



Once again, the ln(A;L) dependence in Eq. ( ]A.28I) does not look like that obtained from 
introducing Eq. f lA.2p into Eq. flA.lip . However the situation here is the same as that 
discussed at the end of Subsection IA.2t leading us to identical conclusions. 



I ^S^CTION A.^ 



Tree-level and one-loop diagrams for Tc^ 



Taking into account the last two Sections, the existence of a perturbative regime and the 
truncation of the series in Eq. f l3.75p at second-order, we can see that just one Feynman-like 
diagram per spectral function of C is necessary to study the tree or loop corrections to these 
spectral functions. In theses cases, just one leading diagram for the one-loop correction to 
Pc^ Fig. IA.2b . as well as one leading diagram for the one-loop correction to -B^ Fig. IA.4b 
are necessary. When applied to T^, this analysis shows that the only diagrams to consider 
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(a) 




Figure A. 5: (a). Tree-level Feynman-like diagram for T^. (b). One-loop Feynman-like 
diagram for T^. The internal dashed lines correspond to two-point correlators of field 
perturbations. 



are the one in Fig. IA.5b for the tree-level terms, and the one in Fig. IA.5b for the loop 
corrections. Such diagrams lead to for T^'^'^'^ and J'^~^°°'p _ 



^f" = iVj^iVj[n^(^2)n^(MP5<^(|k3 + k4|) + 11 permutations] 



27r2 27r2 



27r2 



kl kl |k3 + k4|3 



-|- 11 permutations 



(A.29) 
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+11 permutations 



27r^ 27r^ 27r^ 
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(A.30) 



B 



THE ONE-LOOP INTEGRAL FOR 



We sketch in this appendix the mathematical procedure to estimate the integrals that 
appear when we consider the loop corrections. We only work one integral since the other 
ones are estimated in a similar way. 

The one-loop contribution to the spectrum is: 



47r|k + pPp^ 



-Nl^Vs^ilk + p\)Vs4p) + N;^NI^Vs4\^ + P|)7-,(p) 



1 



(B.l) 



As we can see, the total contribution to ^^"'""p corresponds to three integrals, each one 
having two singularities: one in p = and the other one in p = — k. If the fields spectra 
are scale invariant, the first integral may be written as: 



V. 



1— loop(a) , 



so the actual integral to estimate is: 

/ 



-p2 ivr2 



in 



47r I k + p I '^p^ 



(B.2) 



d^p k^ 



-1 Ik- 



p|3p3 



(B.3) 



This integral is logaritmically divergent at the zeros in the denominator, but there is a 
cutoff at k = L~^. The subscript indicates that the integrand is set equal to zero in 
a sphere of radius around each singularity, and the discussion makes sense only for 
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L ^ <^ k kmax- If we consider the infrared divergences, that means p <C k, we may 
write: ^ 

1= [ ^~47rln(fcL). (B.4) 

To calculate the contribution coming from the other singularity we can make the substitu- 
tion q = k + p. After evaluating this latter integral, we find that the contibution is again 
4:n\n{kL). The integral in Eq. f IB.SP may be finally estimated by adding the contibutions 



of the two singularities: 

d^p k^ 



,^ - 87rln(A;L). (B.5) 
|k + p|V ^ ^ ^ ^ 

More details to evaluate these integrals may be found in Refs. [311 1134[ 1138] . 

The technique to evaluate this kind of integrals when considering vector fields is the same, 
although the procedure is algebraically more tedious. Nevertheless, one can finally arrive 
to the same conclusion. A more detailed discussion about this issue will be found in a 
forthcoming publication |208] . 
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